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Progress by the Georgia Institute of Technology's Laboratory for Radio- 
science and Remote Sensing in developing techniques for passive microwave 
retrieval of water vapor profiles and cloud and precipitation parameters 
using millimeter- and sub-millimeter wavelength channels is reviewed. 
Channels of particular interest are in the troposp h e r ic transmission 
windows at 90, 166, 220, 340 and 410 GHz and centered around the water 
vapor lines at 183 and 325 GHz. Collectively, these channels have potential 
application in high-resolution mapping (e.g. , frczn geosynchronous orbit) , 
remote sensing of cloud and precipitation parameters, and improved 
retrieval of water vapor profiles. 

During the period from July 1, 1993 through December 31, 1993, the 
first wideband millimeter-wave (MMW) and sub-millimeter wave (SMMW) 
atmospheric brightness imagery was obtained using the Millimeter Wave 
Imaging Radiometer (MIR) . The data was recorded during flights of the MIR 
on the NASA ER-2 during CAMEX. A 325-GHz radiometer consisted of a 
sukmil 1 imeter-wave DSB receiver with three IF channels at +/- 1 , 3 , and 8.5 
GHz, and -14 dB DSB noise figure was successfully integrated into the MIR 
for these experiments. 

The submil 1 imeter-wave imagery unambiguously reveals the presence of 
the 325 GHz absorption line as manifested by the spectral signatures 
obtained over clouds and convective precipitation oells, and demonstrates 
the potential for mapping such features using SMMW channels. Analysis of 
the MIR SMMW data along with MIR airborne data fr om cases during TOGA/OQARE 
is in progress, with the primary effort directed toward verifying both 
clear and cloudy radiative transfer models at Y94W and SMMW frequencies. The 
analyses use coincident data from several ground-based, airborne, and 
satellite sensors, including the NASA/MS PC AMFR, the MIT MIS, the EMSP 
SSJVT-2 satellite, collocated radiosondes, ground- and aircraft-based 
radiometers and cloud Ildars, airborne infrared imagers, solar flux probes 
and airborne cloud particle sanpling probes. 

A stannary of ccrparisons between computed and measured clear-air 
brightness tenperatures at the millimeter wave channels 89, 150, 183+/- 
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and 220 GHz has been submit ted for publication. The cl ear - air data 
shew good general agreement between measured brightness temperatures and 
ccnputed brightness tenperatures based on coincident radiosondes, ail though 
the discrepancies between these brightness temperatures show that the 
humidity probes on AIR and VIZ type radiosondes provide readings that are 
too moist in dry regions of the atmosphere. A study using clear-air 325-GHz 
airborne data from CAMEX and ground-based data observed at Georgia Tech is 
underway. 

In order to resolve questions concerning the absolute calibration of 
both the MIR and similar microwave sounders, our study of the scattering 
and emission from microwave blackbody calibration loads has continued. An 
extension of the coupled wave method to two-dimensional periodic structures 
is being developed, and a steady state temperature field solution for the 
wedge- type structure has been developed. The goals of these numerical 
st ud ies are to be able to predict the emitted brightness of such loads when 
used as wideband calibration targets. 

The ability to remotely sense oceanic wind direction using passive 
polar imetric observations is also being studied. A millimeter-wave 
geometrical optics (GO) m od e l for small— amplitude (1— cm peak— to— peak.) water 
waves has been successfully corroborated by 92-GHz laboratory measurements 
of polar imetric e miss ion from a wave tank. A discussion of the GO and 

utility of the third Stokes parameter Ty c ^e<E v E^ 1 *> for passive 
sensing of water wave direction along with a summary of the lab 
measurements has been submitted for publication. Constant bank-angle scans 
of the ocean using a side-looking 92-GHz radiometer on the NASA DC-8 during 
T0GA/O0ARE shew wind-related anisotropies in the upwelling brightness that 
also appear to be of geometrical optics nature in origin. Reduction and 
archival of the 92-GHHz TOGA/OOARE data is nearly ocnplete. 
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Activities within the period from July l, 1993 through VoaeBbote 31, 
1993 under NASA grant NAG 5-1490 have focused on the integr at ion of the 
three-channel 325-GHz receiver on the MIR, operation of the MIR aboard the 
NASA ER-2 during GAMEX (September, 1993 ) , und walifcraUion of the MIR using 
an optimal nonlinear filtering technique, in addition, MIR data analyww 
have continued or commenced, including: (1) several TOGVOCARB 

roeteorological cases, (2) a preliminary analysis of the 325-OJz data, and 
(3) final analysis of dear-air MV brightness data tint were measured 
prior to TOGA/COARE [JacJcscn and Gasiewski, 1993, Appendix Aj . 

The MZR is a joint projec t between the NASA Coddbrd Space flight Center 
and the laboratory for Kadioscience and Remote sensing at Georgia Tech. 
Post Georgia Tech attributions to the MIR end its related scientific uses 

have included basic system design studies, performance analyses, circuit 

and radiometric load des i gn , in-flight software, operation of the MIR 

during field deployments, development of post— fl ight data display arrl 
calibration software, and radiative transfer studies of the. efforts of 
clouds, water vapor, and precipitation on MIR brightness temperatures. 
During Deoaitoer 1993 , th« MTR was transported to Geacg i a Tech for 
characterization and ground-based experiments involving the 32&-<3& 
radiometer. 

In addition to the MIR-related activities, data fl xa a 92-Qfe 
palarinetric radiometer that was deployed an the NASA DC-6 rfm- j pg 
TOGA/COAKE has been calibrated and analyzed. The polarime tr ic capabilities 
of the 92-GHz radla reter were added by Georgia Tech to Investigate the 
relationship between polarized microwave bright n e ss teeperaturee and both 
oriented at on oap n aclc ice dead particles and anisotropic ocean surface 
features. To inpruve the aocuraoy of fully polarimetric radiometers, we 
have designed and begun to construct a prototype di g ital correlator with up 
to 1-G& bandwidth. The correlator will have direct applications in 
airborne polarimetric radiometers suUi at> Us* HAsyMsrc aIvoimU nUawwove 
Precipitation todienster (AMIR) . 
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1. MIR 325- GHz Imagery 

In August 1993 a 325-GHz three-channel DSB radiometer was integrated 
into the MIR. The 325-GHz receiver was built by the ZAX Millimeter Wave 
® ar P ora ^- on according to radicmetric specifications defined by Georgia 
Tech. Die receiver uses a fundamental-mode «i»r with DC bias and a 
frequency-tripled 108 -GHz Gunn oscillator. Three IF channels, one each from 
500-1500, 2000-4000, and 7000-10000 MHz provide spectral sampling similar 
to that of the 183-GHz MIR receiver. All channels are of the total power 
type. With the inclusion of this radiometer the MIR now has nine channels 
at 89, 150, 183+/-1,3,7, 220, and 325+/-l,3,8.5 GHz. 

The measured noise figure of the 325-GHz receiver is -14 dB. Alt hough 
this is relatively high cxmpared to the other four MIR receivers, it yields 
integration noise levels of *2—3 K for a 6 — msec pixel; these have been 
shewn to be acceptable for initial scientific evaluation of the 325.153 GHz 
water vapor absorption line. 

MIR data is being analyzed from several flights during four ER -2 field 
deployments, summarized in Table 1 . In addition to the case studies 
outlined previously 1 , the following new case studies involving the 325-Oiz 
channels have been identified and are being analyzed: 

(1) Clear-air imagery observed during CAMEX (Figures 1 and 2, non- 
cloudy regions) shew brightness temperatures that are, on average, 
approximately the same over all three 325-GHz channels. This is in contrast 
to radiative transfer predictions b ased on the Hebe water vapor absorption 
model. Under the Liebe model [Gasiewski, 1992), the computed brightness 
temperatures for the three 325-GHz channels are nearly identical to those 
of the analogous three 183-GHz channels, and span a range of -25 K. This 
d i screpancy has not yet been explained, but is of critical importance in 
assessing the potential for water vapor sounding at 325 GHz. 

(2) curing several CAMEX flights the MIR observed strong convection 
over both ocean and land. Strip map images of briefness terrperature during 
the flight lines in Figures 1 and 2 show several interesting features of 
the 325 GHz imagery. First, virtually all of the convection that was 


1 Semiannual Status Report #5 f or NASA grant NAG 5-1490, June 30, 
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detected in the channels at 220 GHz and lower in frequency was dsbecfaed 
with comparable brightness variations at 325 (3tz. Thus, water vapor 
screening of convective cell tops does not seen to be significant. This 
observation contrasts with sene previous hypotheses which suggested that 
water vapor opacity wuld preclude detection of cells at 325 (3iz. 

Second, the 325— GHz channels respond to clouds and pr ecip i tati on in a 
mcnotcnic fashion. This response is a resul t of the radionetrically warn 
background caused by the opaque lower atmosphere. In contrast, the 89 and 
150 GHz channels exhibit snail increases in brightness in response to thin 
clouds; the brightness then decreases as cloud opacity and scattering 
increase. In addition, the scattering and absorption at 325 GHz are larger 
than at 89 or 150 GHz. For these two reasons the 325 GHz channels are 
significantly more sensitive to thin clouds and weak precipitation, as can 
be seen by oonparing the apparent edges of rain cells in the 325 and 89/150 
GHz channels. 

Third, the 325-GHz channels unambiguously show the presence of the 
325.153 GHz water vapor absorption line. This is evidenced by the increased 
cloud sensitivity seen in the 325 +/-8.5 GHz channel relative to the 325+/- 
1 GHz channel. Increased absorption by water near the line center reduces 
the reflectivity of the cell top, thereby raising the brightness 
temperature. The same phenomena is seen in the 183+/-1,3, and 7 GHz 
channels. 

Applicable supporting data from other instruments is being ccrpiled for 
both the above CAMEX 325-GHz analyses and ongoing TOGA/CQARE analyses. 
These data are from the Massachusetts Institute of Technology's Millimeter- 
wave Temperature Sounder (MIS) , the NASA/MSPC Advanced Microwave 
Precipitation Radiometer (AMFR), the NASA Airborne Ocean Color Imager 
(Add) , the EMSP SSVT-2 radiometric sounder, the NASA/GSPC Raman water 
vapor lidar, and collocated radiosondes. 

2. MIR Calibration 

In order to aooermodate the relatively large noise levels of the 325- 
GHz channels, as well as to improve the calibration of all the other 
channels, a nonlinear calibration filter was developed. The filter is based 
on: (1) identification and removed of non-stationary features in the 
single-scan gain and offset data, (2) subsequent optimal time-invariant 
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filtering of the residual gain and offset signals, and (3) reconstruction 
of the overall gain and offset waveforms. The ncn-stationary features 
induce jumps and spurious noise. An iterative technique based on the CLEAN 
algorithm is used for their detection. Upon their removal, an opt im al time- 
invariant linear filter (the Wiener filter) is constructed for each 
radiometric channel. The filters are based on the esti m a t ed noise levels 
and autocorrelation time constants for the particular channel. In this 
manner, the unique statistical characteristics of each channel are 
aooonmodated. 

The nonlinear calibration procedure has been shewn to decr ea se the 
contribution of calibration noise in the brightness temperature imagery , 
particularly for the 325-GHz channels. 2 The implementation of the optimal 
calibration algorithm, accommodation of the new 325-GHz channels, and the 
provision of several new data editing features required major modifications 
to the post- flight data analysis software. Most of these have new been 
completed. 

3. Calibration load Analysis 

The emission temperature of a calibration target must be known to 
better than -0.5 K to be useful in calibrating radiometers. To precisely 
predict the emission, it is necessa r y to de termin e both the physical 
temperature and electromagnetic field distribution within the target, which 
is typi call y an array of absorbing wedges or pyramids. The emission 
temperature is the overlap integral of the thermal temperature and a 
normalized function that is proportional to the electromagnetic power loss 
density. This analysis is particularly important for targets used in 
airborne radiometers since these targets are subject to large thermal 
gradients caused by advective cooling. Precise knowledge of the 
electromagnetic distribution in wedge and pyramid arrays is also important 
in the design of microwave absorbers for a variety of purposes, for 
example, anechoic chambers. 

To this end, a steady state numerical solution to the heat equation for 
a wedge-type structure has been developed using a coupled harmonic method. 


2 The basic technique is described by Adelberg and Gasiewski [1993, 
Appendix B]. 
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The prel imin a r y results of the thermal model suggest that temperature 
gradients near the absorbing (and, hence, emitting) tips of the wedges are 
significant enough to warrant consideration in calibration load analysis 
[Gasiewski and Jackson, 1993] . 

4. Polar imetric Micrwave Radicmetiv 

Polar imetric microwave radianetry has been shown to have potential in 
spaoeborne remote sensing of ocean wave direction [Wentz, 1992; Dzura 1992] 
and, possibly, in detection of oriented thunderstorm anvil ice [Evans and 
Vivekanandan, 1990]. We have investigated this both experimentally and 
theoretically using fully polar imetric laboratory measurements at 92 GHz of 
upwelling emission from a fresh-water wave tank [Gasiewski and Kunkee, 
1993b, see Appendix C; Runkee and Gasiewski, 1993]. The measurements, which 
were m a d e using the NASA/GSFC 92-GHz polar imetric radiometer, 3 are well 
corroborated by a geometrical optics model for anisotropic surface 
emission. Both model and measurements show that significant brightness 
variations in the first three Stokes parameters can be produced by only 
moderately striated dielectric surfaces. 

Specifically, the laboratory measurements show a predictable dependence 
of T u =Re<E v Ej 1 *> on the direction of the water wave, with peak-to-peak 
amplitudes of up to 20 K at steep observation angles. Moreover, the Ty 
angular variations are in phase quadrature with similarly strong variations 
exhibited by T v and T^, suggesting that passive remote sensing of surface 
wave direction can be facilitated by polar imetric microwave radicmetry. 
Both the measurements and model calculations consider all four Stokes' 
parameters, althou^i the last of these (T v =Im<E v E^ 1 *>) is very small and not 
expected to be useful for geophysical remote sensing of the troposphere or 
surface. 

The encouraging results of the wave tank experiment prompt ed an 
airborne field experiment during TCGA/OQARE to investigate passive 
polar imetric measurement of ocean wave direction. Indirectly, one can be 
expected to infer atmospheric wind direction from such measurements. 
Electronic and hardware modifications necessary to operate the polar imetric 


3 The polar imetric capability was added by Georgia Tech under a 
previous NASA grant (NAG 8-829) [see Gasiewski and Kurikee, 1993a] . 
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radiometer on the NASA DO-8 aircraft were performed. A total of twenty DC-8 
flights occurred during TDGA/CQARE, including seven low-altitude (1.5-4 km) 
constant bank-angle maneuvers designed to provide views of the ocean 
surface at a constant observation angle and over a range of azimuthal 
angles [Kunkee and Gasiewski, 1994]. 

Post-mission calibration of nearly all of the 92-GHz polar imetric data 
has been completed. A plot of the constant bank-angle data for an incidence 
angle of 65° (Figure 3) shews residual peak-to-peak brightness variations 
over azimuthal angle of amplitude -3 K for T v and T^. The shape and 
amplitude of the variations are closely related to those found by Wentz 
[1992]. Even for observation at nadir (Figure 4), sinusoidal brightness 
variations of amplitude -1 K (complementary in T v and T^) can be seen. 
Although the nadir variations are not large enough for retrieved purposes, 
they are large enough to adversely impact brightness measurements for 
atmospheric sounding purposes, for example, wet path delay measurements or 
water vapor sounding. The TDGA/CQARE data, along with the laboratory wave 
tank m e a surements corroborate the hypothesis that anisotropic ocean surface 
signatures can be observed using microwave radiometers. 

In addition to striated water surfaces, it has been hypothesized that 
polar imetric microwave signatures in T v , T^, and T\j will be produced by 
oriented ice particles, for example, in electrified cirrus anvils. Indeed, 
microwave depolarization signatures in spaoe-to-greund communications links 
have been associated with li^itning discharges [Cox and Arnold, 1979]. 
Although polarized microwave signatures from hydremeteors have not yet been 
definitively observed in the TDGA/CQARE data, we sure analyzing the origin 
of apparent T v -T^ differences of up to +/-10 K observed at 92 GHz over 
oceanic convection. Supporting data from the other side-looking radiometric 
instruments that flew on the DC-8 is being sought to verify whether or not 
these differences are instrumental. 

5. Digital Correlator for Polar imetric Radiometrv 

In anticipation of the need for precision airborne and spacebome 
polar imetric radiemetry, a high-speed digital correlator for proposed use 
in the NASA/MSPC Advanced Microwave Precipitation Radi ome ter (AMFR) is 
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being developed [Gasiewski and Kunkee, 1992]. 4 A p ro totype A/D converter 
operating at 800 Ms/sec and using standard emitter-ooupled logic (BCL) has 
been demonstrated (Figure 5) . Two such A/D converters along with three high 
speed digital counters are the essential components of the correlator. The 
reason for using digital correlation is that calibration of the cross- 
oorrelaticn channel can be accurately performed using only the standard hot 
and cold calibration targets found in ccnvecticnal radiometers. The 
performance of our pro t otype converter suggests that the n ecessarily wide 
IF bandwidths required for Earth remote sensing (in this case, up to 400 
MHz) can be obtained using the digital technique. Measurements of converter 
hysteresis and development of the ECL counters (to be used in conjunction 
with the A/D converters to form the complete correlator) are in progress. 


4 A proposal for ouch more extensive development of digital 
correlation polar imetry and its evaluation on the NASA/MS FC AMFR has been 
submitted by this PI to NASA Headquarters under the title "Passive 
Measurement and Interpretation of Polarized Microwave Brightness 
Temperatures, " September 1992, NASA control # 2916-RD-074. 
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SMRHY AND FUNS FOR FUTURE WORK 


me investigations conducted over the past six norths have served 
Purposes, including: (l) the first dencnstretim of radirsetric ijnagin, for 
meteorological purposes at 325 Oiz, (2) the development and iaplenentati® 
of cptraal nonlinear filters for radiometer calibration, ( 3 ) identification 
of several case studies for radiative transfer modelling at both 
millimeter- and subnillimeter frequencies, and in both clear and cloudy 
atmospheres, and (4) demonstration of brightness temperature anisotropies 
caused by ocean saves at 92 GHz. future plans shich build on the timings 
of these investigations are outlined below. 

■L. MIR Observations an d Data Analyse 

An ensemble of clear-air 325-GHz observations for the purpose of 
radiative transfer model validation will be cerpiled from the CAMEX data. 
Ground based data measured with the MIR at Georgia Otech will be used to 
supplement this. These data will be corpared with computed brightness 
tenperatures based on collocated radiosondes and the Liebe water vapor 
absorption model [Liebe, 1985]. Steps to inprove the absolute calibration 
of the MIR (see Section 2 below) will be necessary in these cotparisons. 

Several convective case studies fran CAMEX are also being cerpiled for 
statistical analysis. Curing the convective overflight on October 5 , 1993 , 
the ER-2 Doppler radar (EDOP) obtained coincident nadir reflertivit^ 
profiles. This is the first joint EDOP and MIR data set. Because of the 
coincident EDOP radar truth, analysis of the MIR data fran this flight is 
of great interest. 

Two of the MIR case studies identified during TOGVOOARE are 
particularly interesting for radiative transfer experiments, namel y, the 
overflights of cyclone Oliver and the Kavieng ground observation site. To 
this end, we plan to simulate ocean surface and atmospheric conditions 
within the eye of cyclone Oliver to determine the consistency of wideband 
radiative transfer models using the observed brightness data. The radiative 
transfer model will consider the effects of increasing humidity, 
precipitation, and ocean roughness near the eyewall to determine the 
relative contributions of these oorponents to the warm ring. This 
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experiment will require sane data frcm the NASA/MSPC AMFR and JPL APMAR, 
both of which are available and are being requested.** 

Bie Kavieng overflights will be used to provide data for simulating the 
effects of cl cuds on upwelling M»W brightness temperatures. Ground-based 
radiometer data will be used to determine total water vapor ard cloud 
content, while ground-based lidar will be used to determine cloud bottom 
altitudes. R ad iosondes will be further used to constrain the ver tical 
distribution of water vapor. Of interest is a comparison of oomputed and 
observed upwelling brightness temperature, and in particular, the inpact 
that clouds have on these temperatures. 

MIR data is currently stored in raw format on high-density 8 -cm tapes. 
To facilitate meteorological data analysis and to provide a practical means 
of disse m i n ating MIR data to collaborating investigators, the display 
software is being modified for disk-based storage. This software, which 
will be available for use on PC's, will also incorporate a variety of 
simple interactive features for MIR data analysis. The MIR data is aisn 
being archived on 35-nm slides for graphical storage and dissemination. 

2. Radiometer Calibration 

In order to better characterize the RF response and absolute accuracy 
of the MIR, we plan to conduct several tests, including: (1) RF passband 
response measurements using a plasma discharge noise source and IF spectrum 
analyzer, (2) local oscillator interference and reflection measur ements 
using a stepped reflecting plate, and (3) calibration-load foam reflection 
and transmission m e as urements. These simple measurements will provide 
answers to questions concerning the calibration of the MIR and the use of 
the 183 and 325-GHz data in radiative transfer interoenparisons. Particular 
attention is being paid to the 150 GHz channel, which consistently returns 
brightness temperatures that are 5-10 K oolder than expected from an 
absorber im m ersed in liquid nitrogen. 

We are new beginning to use the nonlinear calibration technique for 
operational calibration of all MIR data, including flights during CAMEX and 
TOGA/ CHARE. We plan to archive the calibrated data so that it will be 


5 Collaborating on the TOGA/OQARE data analysis is Dr. J. Vivekanandan 
of the university of Colorado. 
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available to investigators collaborating cn TOGyOARE studies. A 1-GB hard 
disk drive has been obtained under this grant to suppor t the calibration 
and archival effort. Further improvement in the nonlinear calibration 
algorithm is expected to be found by using a lmitiHiTn^n g-innai wiener 
filter. This will be tested, and if successful, will be implemented prior 
to archival. 

Accurate absolute calibration of the MIR (and s imilar radiometers) 
requires that the total reflectivity of the hot and cold loads be less than 
”1%, and known to better them 0.1%. Manufacturer's specifications typically 
provide only the specular ccnponent of the reflectivity , which is thought 
to be substantially less than the total reflectivity. In order to refine 
the MIR calibration, we plan to extend the study of the electromagnetic 
characteristics of wedge-type blackbody loads to the more desirable 
pyramidal loads. 

Cur approach is to develop numerical models for one- and two- 
dimensional ly periodic lossy gratings using the coupled wave method. We 
currently have software based on the coupled wave method to predict the 
reflectivity of one-dimensional ly periodic loads of arbitrary dielectric 
profile. Extension of the coupled wave model to two-dimensianally periodic 
surfaces appears feasible, and will be continued. In addition to the 
electromagnetic analysis, work on the steady state thermal analysis for 
two-dimensionally periodic calibration loads will be continued. Although, 
the radicme trie calibration load reflectivity analysis is of importance in 
underst andin g precision radiometer calibration, it is of secondary 
importance relative to the M4W and SM4W data analysis. 

3. Pol arime trie Radicmetry 

Further investigations of polar imetric radicmetry will consist of 
follcw-up analyses of 92-GHz polar imetric data observed during TOGA/OQARE, 
and the development of a statistical ooean surface emission model using 
geometrical optics theory for corroboration of the measured TOGft/OQARE 
data. The model will includ e the effects of ocean surface foam. The 
limitations of the gecmetrical optics model in explaining the wave tank and 
ocean surface data will be of interest; the validity of this model at 92- 
GHz will provide insight into the applicability of geometrical optics-based 
models for higher microwave frequencies. Information on the polarizing 
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properties of surfaces will be useful in remote sensing of ocean surface 
characteristics and in understanding the inpact of ocean surface emission 
on passive atmospheric sounding. 

To reduce the oonplexity of calibrating a polar imetric radiometer, the 
development of a working digital cross-correlator will be continued. The 
cross correlator will be a protot y pe for the one proposed to be used on the 
NASA/MSFC AMFR. The bandwidth of the cross-correlator will be approximately 
500 MHz, making the device useful for wideband radiometric channel studies. 
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Table 1. 


MIR integration and data flights: 1992-93. 

Sortie # Date Time Instruments 

(UPC) & Remarks 


Moffett Field, CA: 

92-087 5/11/92 1900-2200 

92-089 5/14/92 2315-0515 

92-090 5/15/92 2320-0500 

Wallops Island, VA: 

92-130 7/23/92 2100-2315 

92-131 7/29/92 0700-1330 

92-132 7/30/92 0700-1330 

92-134 8/2/92 0700-1030 

92-135 8/3/92 0700-0615 

92- 140 8/6/92 0700-1330 

TOGA/OOARE (Townsville, ADS) 

93- 053 1/12/93 2130-0430 

93-054 3/17/93 2300-0710 

93-055 1/19/93 0130-0922 

93-056 1/25/93 2315-0700 

93-057 1/30/93 

93-058 1/31/93 2200-0600 

93-060 2/5/93 1430-2050 

93-061 2/7/93 1555-2115 

93-062 2/9/93 1815-0025 

93-063 2/10/93 1430-2225 

93-064 2/19/93 

93-065 2/20/93 1900-0335 

93-066 2/22/93 1900-0205 

93-067 2/24/93 2000-0315 

CAMEX (Wallops Island, VA) : 
93-164 9/12/93 

93-165 9/15/93 

93-166 9/19/93 

93-167 9/25/93 1751-2002 

93-168 9/26/93 1908-2341 

93-169 9/29/93 0106-0458 

93- 178 9/30/93 2012-0213 

94- 001 10/3/93 

94-OOx 10/5/93 


MIR, MTS (u) ,AOCI 
Mm,MPS(d),AOCI * 

MIR, MIS (U) ,AOd * 

MPS(d) 

MPS(d) +* 

MTS(d) +* 

MPS(d) + 

MTS(u) + 

MIS(u) +* 

Radiation (93-01-06) 
Convection (93-01-07) 
Convection (93-01-08) 
Radiaticn/Kavieng ( 93 - 01 - 09 ) 
Pilot proficiency 
Radiation (93-01-10) 
Convection, Oliver overflight 
(93-01-11) * 

Cyclone Oliver, MPS(d) 

Oliver overflight (93-01-13) 
Kavieng, MES(d) (93-01-14) * 
Pilot proficiency, MIS(u) 
Convection (93-01-16) 
Convection (93-01-17) 
Radiation (93-01-08) 

Convection, Transit Flight 
Flight #1 Eng. Test (aborted) 
Flight #2 Eng. Test 
Eng. Test Flight 
(times unsync'd) 

Flight #3 
Flight #4 (AIRS) 

FLight #5 (SStyT-2) 

Flight #6 (convection) 

Flight #7 (convection, EDOP) 
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Notes for Table 1 


* SSH/T -2 satellite underpass. 

+ Ground-based Raman HoOy lidar overflight, 
(xx-xx-xx) indicate coincident DC-8 flights. 

Note: "u" or "d" indicate i$>-looking or down- 
looking, respectively. 

All CftMEX flights include the 325 -<3iz channels. 
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Figure 1 (p. 19) : Strip map images of oceanic convection observed by 
the MIR on the ER-2 at 20-km altitude (CAMEX, 9/26/93, 2035-2053 
UTC) . Several snail convective cells indicated by brightness 
temperatures lower than 150 K are shown. 


Figure 2 (p. 20) : Strip map images of oceanic convection observed by 
the MIR on the ER-2 at 20-km altitude (CAMEX, 9/26/93, 2058-2116 
UTC) . The anvil region of the raincell near the cursor (2104) is 
clearly outlined by the 325 GHz channels. 
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Figure 3. Vertical (A) and horizontal (B) brightness variations at 92 
GHz for constant bank-angle flight over ocean during HXyyOQRRE (2/8- 
9/93). The altitude is 1.5 km, and the radiometer beam angle is 65° 
from nadir. The abscissa is the angle of the beam relative to the 
estimated direction of the ocean waves. 





Wave Azimuthal Angle (degrees) 


Figure 4. Vertical (solid) and horizontal (dashed) brightness 
variations at 92 GHz for constant bank-angle flight over ocean during 
an integration flight in preparation for TOGA/CQARE (12/12/92). Ihe 
altitude is 1.3 km, and the radiometer bean angle is approxi ma tel y at 
nadir. Ihe abscissa is the angle of the beam relative to the 
estimated direction of the ocean waves. Ihe vertical and horizontal 
brightness temperatures are antiphased and of second harmonic 
variation. 
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Abstract 


A comparison of clear-air brightness temperatures is performed between radiometric mea- 
surements and atmospheric radiative transfer calculations. The measurements were made 
using the NASA Goddard Space Flight Center’s Millimeter-wave Imaging Radiometer 
(MIR) in a series of airborne and ground-based atmospheric experiments at six millimeter- 
wave frequencies: 89, 150, 183.3±1,±3,±7, and 220 GHz. With the inclusion of the 220 GHz 
channel, these measurements are the first passive observations of the atmosphere made 
simultaneously at the six frequencies. The MIR was operated concurrently with supporting 
meteorological instruments (radiosonde and Raman lidar) to construct a paired set of both 
spatially and temporally coincident calibrated brightness temperatures and atmospheric 
profile parameters. Calculated brightness temperatures based on the measured atmospheric 
profile parameters were obtained using a numerical radiative transfer model. Incremental 
water-vapor weighting functions were used to study the impact of radiosonde hygrometer 
errors on the radiative transfer calculations. The aircraft-based brightness temperature 
comparisons are generally good for the channels sensitive to the lower atmospheric levels 
(89, 150, 183.3±7 and 220 GHz), but show discrepancies of up to 11 K for the opaque chan- 
nels (183.3±1, and ±3 GHz) caused primarily by radiosonde bias. The ground-based calcula- 
tions are similarly found to be sensitive to hygrometer errors in the lower atmosphere. 
Ground-based comparisons between MIR observations and lidar-based calculations are typi- 
cally within ±6 K. 
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1.0 Introduction 

In constructing water-vapor profile retrieval algorithms using spacebome passive 
millimeter-wave observations the relationship between the profiles and upwelling radiation 
must be known. Radiometric ch ann els for water-vapor profile retrieval include those near 
the weak 22.235 GHz water-vapor absorption line and the stronger 183.310 GHz line. 
Recently, the potential uses of the 325. 153 GHz water-vapor line have also been identified 
for the spaceborne retrieval of tropospheric water-vapor profiles and cloud parameters [1]. 
In developing retrieval algorithms using these fines, verification data i6 required consisting 
of radiosonde measurements of temperature, pressure, and water-vapor density profiles. 
However, the accuracy of a candidate water-vapor retrieval algorithm is limited by the accu- 
racy of the radiosonde data. Comparisons between clear-air multispectral brightness tem- 
peratures from a radiometer and radiative transfer calculations can be used to investigate 
the impact of radiosonde errors on retrievals. 

Currently, tropospheric sounding facilities in the continental United States use at 
least three different radiosonde packages: AIR, VIZ, and Vaisala. Both radiosonde hygrome- 
ter errors [2] and differences in relative humidity reporting practices [3] have been noted. 
Such errors at specific altitude levels will influence brightness calculations for channels sen- 
sitive to that level through the relative humidity weighting functions. These functions 
express the response of a radiometric observing system to small perturbations in the water- 
vapor density profile. For example, for spaceborne clear-air observations within 3 GHz of 
the 183.310 GHz water-vapor line, the weighting functions peak at altitudes of 7 km and 
above, where the corresponding relative humidity values range from 0 - 15%. Because a 
radiosonde hygrometer’s sensitivity is monotonically related to the number density of water 
molecules present, their accuracy degrades in this region of the troposphere. Brightness 
comparisons at these channels can be expected to show this degradation as an increased dis- 
crepancy. 

Clear and cloudy brightness temperature comparisons between radiometric mea- 
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surements and radiosonde-based radiative transfer calculations were performed by Westwa- 
ter et al. [4] and England et al. [5] for ground-based vertical-angle observations and for 
frequencies up to 90 GHz. Mid-altitude ( < 9 km) aircraft-based comparisons were performed 
at 89 and 157 GHz by Foot et al. [6] (zenith- and nadir-directed) and English et al. [7] 
(zenith-directed) to study water-vapor continuum absorption. For these frequencies the 
water-vapor weighting functions show no appreciable sensitivity to water-vapor perturba- 
tions in the upper troposphere. Detailed satellite-based comparisons at 89, 150, and 
183.3±1,±3,±7 GHz have also been performed by Falcone et al. [8] and Morrissey et al. [9], 
but were impeded by the sensor’s relatively large field-of-view (-50-80 km) and the subse- 
quent difficulties in modeling the observed background and the horizontal variability of 
water-vapor density in the atmosphere. In the above investigations, the effects of radio- 
sonde “dry-end” upper atmosphere reporting errors on the theoretical calculations were not 
reported. 

In this paper, we compare the observations made with the NASA Goddard Space 
Flight Center’s Millimeter-wave Imaging Radiometer (MIR) and radiative transfer calcula- 
tions based on simultaneous and collocated radiosonde and Raman lidar measurements. 
The MIR is a cross-track scanning aircraft-based radiometer currently operating with the 
following six millimeter-wave (MMW) channels: 89, 150, 183.3±1,±3,±7, and 220 GHz, and 
with provisions for three submillimeter-wave channels at 325.2±1,±3, and ±8.5 GHz. In this 
investigation the MIR flew aboard NASA’s high altitude (-20 km) ER-2 aircraft during 
deployments out of the NASA Wallops Flight Facility, VA, and the NASA Ames Research 
Center, MofTet Field, CA. Brightness comparisons are made over both land and ocean back- 
grounds. Water-vapor weighting functions are developed for the MIR’s frequencies and for 
both zenith- and nadir-directed observation geometries. Results of ground-based zenith- 
directed comparisons are also discussed. A preliminary study has been done by Wang et al. 
[10] for a subset of the above aircraft-based MIR comparisons. 
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2.0 Millimeter-wave Radiative Transfer and Radiometry 


2.1 Radiative Transfer Theory 

The absorption and emission of microwave radiation by the atmosphere exhibits res- 
onances due to transitions among the quantum energy states of several molecular constitu- 
ents. The water-vapor absorption resonance at 183.310 GHz, and the oxygen resonance at 
118.750 GHz, are most important in the frequency range of this study, as shown in Fig. 1. 
The absorption near these opaque lines, as well as the continuum absorption near the win- 
dow frequencies of 90, 150, and 220 GHz, is a function of the local pressure (P), temperature 
( T ), and the density (p) of the absorbing constituents. Absorption models (e.g. [11], [12]) 
relate vertical distributions of P, T, and the constituent p’s to profiles of the absorption coef- 
ficient, a(f,z ) , where f is the frequency and z is the altitude. The integrated radiative 
transfer equation (IRTE) relates the observed radiant energy to profiles of absorption coeffi- 
cient and the intensity along a propagation path to the observer [13]. The radiation field is 
typically characterized by a brightness temperature, defined by scaling the radiation inten- 
sity (in W/m 2 -St-Hz) by A /2k, where A is the wavelength and k is Boltzmann’s constant. 
Using this definition, the brightness temperature for a blackbody radiator at physical tem- 
perature T is: 


T B lf) = 


hf 


k(e hf/kT -l) 


(1) 


where h is Planck’s constant. Note that this definition for brightness temperature differs 
from the Rayleigh-Jeans approximation [14], wherein Eq. (1) is expanded to first order in 
hf/kT and leads to Tg (f) - T. 

For a clear and horizontally-stratified atmosphere, the brightness temperature 
observed at angle 6 relative to the zenith for downward-looking platforms at height h is: 
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h 

T B (f,h,6) = jT Bt (f,z) a(f,z)sec6 e~ T ^' th),tc9 dz + e~ x ^’ 0,h)tBce T B0 (f,6) (2) 

o 

for 14 < ( n/2 ) , where 


r B0 (/;e) »(i-r(©))r B< 


(3) 


( - 


\ 


+ r(6 ) 



z)a(f,z)sec6 




V 0 ' 

Tjj, and T CB are the surface and cosmic background brightness temperatures, respectively. 
T Bz (f, z) is the local brightness temperature determined using Eq. (1). The integrated opacity 
over the path interval is: 


h 2 

t (f,h v h 2 ) = j a(f,z)dz (4) 

In Eq. (3), the specular surface reflectivity rid) represents the fraction of power reflected by 
the surface that is copolarized with the observing instrument's antenna. Over an ocean 
background rid) is a function of frequency and surface roughness, whereas over land rid) is 
typically small ( ~ 10%), broadband, and largely independent of angle. For ground-based, 
upward-looking observations, the brightness temperature is: 




(5) 


for |6| > (n/2) . 

Note that the IRTE (Eqs. (2) or (5)) is exact in as far as no approximations are made 
in Eq. (1) for brightness temperatures. However, in obtaining T B (f, h, 6 ) , it is more conve- 
nient to work with expressions involving physical temperatures rather than brightness tem- 
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peratures. Tb this end, we substitute Eq. (1), expanded to second-order in hf/kT , for 
T Bz (f,z) in Eqs. (2), (3), and (5), and for T Rg in Eq. (3). Equation (1) with no approximation 
is substituted for T CB , where 7 = 7 c and T c is the cosmic background temperature 2.73 K. 
As a consequence of these substitutions, T C b is now defined as [15]: 


7 

J CB 2k 


( e hf/kT ‘+l' 
e hf/kT ‘-V 


( 6 ) 


and T B (f, h, 6) now implicitly contains a second-order correction, i.e., in Eqs. (2) and (5) 
T b {f, h, 6) is now interpreted as T B (/“, h, 6) + hf/ 2k. Due to the nonlinearity in brightness 
temperature with physical temperature for temperatures near 7 C , the second-order correc- 
tion is adjusted by T CB from Eq. (6). In using Eq. (1) to second order, T R (/) is accurate to 
-0.15 K for frequencies up to 340 GHz and for atmospheric temperatures as low as 150 K. 
In contrast, the maximum error incurred using the Rayleigh-Jeans approximation over this 

range is 8 K. 


2.2 Incremental Weighting Functions 

The IRTE is nonlinear in the relevant atmospheric parameters (i.e., P, 7, and partic- 
ularly the constituent p’s) due to their influence on the absorption coefficient a, which enters 
the equations in a nonlinear manner. However, weighting functions for particular atmo- 
spheric parameters can be determined by linearizing the IRTE, where the respective param- 
eter is incrementally changed from a nominal value. Incremental weighting functions 
express the brightness temperature response of an observing system to small perturbations 
in the atmospheric parameters of interest. In the present study, we are interested in the 
effects of perturbations in the water-vapor density profile on observed brightness tempera- 
tures caused by radiosonde errors. Incremental water-vapor weighting functions for the spe- 
cial case of vertical viewing angles in downward-looking and upward-looking observing 
systems have been discussed by Schaerer and Wilheit [16], and Westwater et. al. [4], respec- 
tively. 
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Tb linearize the IRTE, we seek brightness changes ST B caused by deviations in both 
the temperature 5T(z) and absorption coefficient $a(z) from their nominal profiles. For 
microwave absorption in the troposphere, the absorption coefficient a(f,z) is a function of 
7Xz) and p„(z), i.e., a = a(T,p v ) . Using Eq. (4), the absorption perturbation leads to an 
opacity perturbation, given to first order by: 


“ (! w+ a7>) & (7) 

where Sz is the thickness of the perturbed layer. The corresponding change in brightness 
temperature to perturbations at height z is: 


ST b (z) = 


* T B 

BT(z) 


ST(z) + 


dT B 

dT(z) 


Sr (z) 


( 8 ) 


From perturbation analysis, the two terms on the right-hand side of Eq. (8) can be deter- 
mined directly from Eqs. (2) or (5), where we use physical temperatures for brightness tem- 
peratures. For downward-looking observations we can write ST B as: 

h h 

ST b = J Wj, (z, 6) ST (z ) sec 6 dz + j ( z , 6) 8p v (z) sec 6 dz (9) 

0 0 

where the incremental temperature weighting function, W T (z, 6 ) , is given by: 


W t (z,6) = a(f,z)eecd ( e -T( ^’ ,ec e + r ( 0) e~ x ^' °’ A) * ec e e ~ x ^' °’ “ c 

+ sec© ^T(z)e~ 1{f,t ' h) “ c ®- ^T(z')a{f,z') sec© e ~ T{f ' g ‘ ,h) uece dz' 


( 10 ) 


- R(z, 6) e-^ 0 *^‘ ece } 
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and the background term B(z, 6) is: 


K..O - (WO)Vr(l|{ir ffl 


(ID 


jT (Z') a <f,z') sec 6 e x{f 0 ’ t ' ) ™ c9 dz' + jr (z') a(f,z') aec6 


-f (f, 0,*') IK 9 


dz' 


} 


The downward-looking incremental water-vapor density weighting function, W p> (z, 6 ) , is: 


W. (z, 6) = sec 
p ’ *P„ 


t 

e ^ r(z)e -T( ^’ z ' A) * ec0 - jVfz') a(f,z') seed e - T V'*' h)tec9 dz' (12) 


- B(z, 6) e 


-r(f, 0, A)*ecS 


} 


For the ground-based upward-looking observations, the incremental temperature 
weighting function Wy,(z, 6) is: 


W T (z, 6) - a (f, z) sec 6 e 


0) »ec 0 


( 13 ) 


da a 
+ — seef) e 

dT 


- Jr(z')a(/,z - )sece e'’ 11 '-” ‘' c "dz } 


and the upward-looking incremental water-vapor density weighting function, ( z , 6 ) , is: 


W n (z, 0) = sec0 e 
p * d P v 


r r(/.z. 0 ) mc T(z)-T cb e - r(f '-' x) ,ecS 

- Jt(z’) a(f,z') sec 6 e ~ r(f ' z '' z) 8ec0 dz']- 


( 14 ) 
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with T C b for both observing regimes determined from Eq. (6). 


The first term in the incremental weighting function W r (z, 6) represents the linear 

response to temperature deviations. The terms proportional to Ba/BT and 3 aJBp v inEqs. 

(10), (13), and Eqs. (12), (14), respectively, describe the IRTE’s departure from strict linear- 
ity with the absorbing constituent in question. These terms also Bhow that the difference in 
functional form of the response to T and p„ is due only to the multiplicative factors Ba/BT 
or Ba/Bp v , respectively. Since a is relatively insensitive to T, Ba/BT is often considered 
negligible (e.g., [17]) and the IRTE is considered linear in T. 

In Figs. 2 and 3, incremental weighting functions ( z , 6) for aircraft- and ground- 

based geometries are shown for U.S. standard summer mid-latitude atmospheric conditions 
[18] using an exponential water- vapor density profile with a 2-km scale-height and a surface 
relative humidity (SRH) of 75%. The plots are of relative humidity weighting functions, 
which are related to water-vapor density weighting functions by: 

w s <2.e, - w p .( 2 ,e)^ as) 


where p M( (z) is the saturation water-vapor density [19]. Figure 2 shows downward-viewing 
(at nadir) relative humidity weighting functions computed over both an ocean (Fig. 2a) and 
land background (Fig. 2b) at several MMW frequencies. In Fig. 3 ground-based zenith-view- 
ing weighting functions are shown for the same frequencies. Given an assumed relative 
humidity perturbation SR over thickness &, brightness temperature perturbations can be 
determined directly from these figures. 

As expected, the peaks of the 183-GHz downward-viewing relative humidity weight- 
ing functions move lower into the atmosphere for frequencies progressively farther from 
183.310 GHz resonance. For example, observations with the 183. 3±1 and 183. 3±7 GHz 
channels are most sensitive to radiation originating from altitudes of 8-12 km and 3-6 km, 
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respectively. The sign of window-channel weighting functions can change, as seen in their 
response over a land background (Fig. 2b) and in the 220 GHz channel over ocean (Fig. 2a). 
Therefore, positive-valued hygrometer errors in the atmospheric levels above 3 km can 
decrease brightness temperature calculations, and in the lower levels can increase bright- 
ness temperatures. The weighting functions are dependent on the unperturbed nominal 
temperature and relative humidity profiles. However, if the SRH deviates from 75% (consis- 
tent with the radiosonde profiles collected for this study) by ±20% RH the magnitude of the 
weighting function peaks and their positions are not significantly changed. Note that for 
upward-viewing observations (Fig. 3), the 150 GHz channel is the most sensitive to water- 
vapor perturbations in the lower atmosphere. 

2.3 Millimeter- wave Radiometric Principles 

For an incident unpolarized brightness field T B (f, h, £2 ) , where £2 is the propagation 
direction, the antenna temperature T A of a radiometer can be expressed in terms of the radi- 
ometer’s antenna gain pattern GU7) and normalized receiver frequency response Hif) as: 

•o 

r A (h)~ iff (f> i 2 1 a (fl> t b if, h, a) dadf as) 

0 4 * 

where fi = J~| \H{f)\ 2 df is a spectral normalization factor and h is the radiometer’s height. 
Thus, the measured temperature T A is a weighted average of the brightness temperature 
reaching the receiver from all directions in space and from all frequencies in the receiver 
passband. In practice, the bulk of the received energy for most radiometers enters through 
the relatively narrow main beam of the antenna. 

The underlying atmospheric thermal emission observed by the radiometer is a ran- 
dom process. The resulting measurements are estimates of the power in this process and are 
themselves a random process. Consequently, the radiometer output fluctuations are propor- 
tional to the sum of T A and the receiver noise temperature, T^. The resulting sensitivity, or 
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minimum detectable temperature, of a total-power radiometer is given by: 



where B is the IF channel (3-dB) bandwidth, ris the post-detection integration time, and 
( Og/(g )) 2 is the relative predetection gain variance. The gain variations occur at time 
intervals much less than a typical calibration interval (-1-2 sec). Tb minimize the error in 
the estimate of T A , it is desirable to increase the radiometer bandwidth (£) and integration 
time (t). However, for practical measurements, the bandwidth is limited by the need for both 
detailed spectral information and interference rejection, and the integration time is limited 
by the required temporal sampling time. 

2.4 Millimeter-wave Imaging Radiometer Description 

The Millimeter-wave Imaging Radiometer (MIR) is a total power six -channel imag- 
ing spectrometer designed for airborne studies of passive microwave retrieval of tropo- 
spheric water vapor, clouds, and precipitation parameters [20). The MIR is a cross-track 
scanning radiometer and can be configured for either airborne nadir-viewing aboard the 
NASA ER-2 high-altitude aircraft or ground-based zenith-viewing. Table I summarizes the 
MIR channel design specifications. The radiometric sensitivities were calculated using Eq. 
(17) for a 75-ms integration period, relative gain variance of lO - ®, and nominal antenna tem- 
perature of 200 K. 

The MIR consists of a receiver housing and data acquisition system. For aircraft- 
based observations, the receiver housing mounts into the camera port of the ER-2 superpod 
nose cone. Within the receiver housing an ambient pressure cavity contains the scanning 
mirror, stepper motor, hot and cold calibration targets, and temperature sensors. A pressur- 
ized enclosure houses the lenses, feedhoms, mixers, IF and video hardware, and data acqui- 
sition analog and digital electronics. A PC-AT compatible computer system, mounted above 


11 



the receiver housing, executes configuration and control programs stored in a nonvolatile 
hard disk emulator. Digitized radiometric data received from A/D boards and system opera- 
tional and diagnostic data are recorded on 8mm EXABYTE magnetic tape cartridges. 

Each of the five MIR feedhoms is located at the focal point of a corresponding lens 
integrated into the pressurized receiver housing wall. All lens-feedhorn pairs have nearly 
identical 3.5° beam widths. The five mixers are operated in a double-sideband configuration: 
the 89 and 150 GHz mixers are balanced; the 183.3, and 220 GHz mixers are sub-harmoni- 
cally pumped. All LO’s are InP Gunn oscillators where the second harmonic is extracted and 
used in the mixing process. After IF amplification, the radiometric signals are bandpass fil- 
tered. For the 183.3 GHz mixer, a triplexer is used to divide the IF signals into three chan- 
nels. Detected signals are amplified, low-pass filtered, and subsequently, sampled by 12-bit 
A/D converters. DC feedback is added to each video signal via an 8-bit D/A converter. This 
permits in-flight adjustments to compensate for slowly varying drifts in receiver gain and 
noise temperature. The feedback can compensate for receiver gain variations of ±2 dB and 
receiver noise temperature variations of ±20%. The feedhoms are directed, in turn, to hot 
(-325 K) and cold (-240 K, ambient temperature) calibration targets that serve as known 
blackbody radiometric sources. Absolute calibration accuracy is estimated at better than 1 K 
from knowledge of both the temperature sensor accuracy (better than 0.1 K) and the electro- 
magnetic scattering characteristics of the targets [21]. In-flight observations of both calibra- 
tion targets indicate channel sensitivities consistent with those shown in Table I (Eq. (17)). 


3.0 Description of Experiments 

3.1 Wallops Experiments, July-August, 1992 

Simultaneous and collocated measurements involving the MIR, a Raman lidar, and 
radiosondes were performed on five days between 29 July and 6 August 1992, at the NASA 
Wallops Island Flight Facility, Virginia. Motivation for these experiments included using 
MIR data to validate the calibration accuracy of the Special Sensor Microwave (SSM/T-2) 
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water vapor sounder [8], collecting data for future water-vapor retrieval comparisons, and 
performing clear-air brightness temperature comparisons. Statistics for MIR-observed 
brightness temperatures, over ocean and land backgrounds, are shown in Table II. The sta- 
tistics are computed from successive nadir samples taken over 1-2 minutes during selected 
periods of each flight. For the comparisons, MIR observations and calculations derived from 
radiosonde and ground-based Raman lidar data are used. Both ground-based and aircraft- 
based experiments were made with the MIR. 

Ground-based observations were performed over a 100° swath centered at zenith 
from 2200 to 0300 UTC on 29-30 and 30-31 July 1992. The observations were made during 
clear and stable atmospheric conditions and within 5 km of an active radiosonde network. 
Aircraft-based observations were made from the NASA high-altitude ER-2 aircraft during 
overflights of both coastal regions near the Wallops Island facility, and ocean flight-tracks, 
parallel to the coastline, approximately 200 km offshore. Three 6.5-hour and two 3-hour 
flights were made during local nighttime conditions, as summarized in Table III. A total of 
27 Wallops radiosonde profiles were used in the clear-air comparisons. Temporal and spatial 
collocation differences were typically better than 45 minutes and 20 km, respectively, 
between ER-2 positions, radiosonde launches, and Raman lidar observations. Data from the 
Cape Cod flight segments were used for SSM/T-2 calibration studies. 

A typical brightness temperature time-series is shown in Fig. 4, for the 29 July flight 
from 0900 to 0917 UTC (ocean to land flight track), and 0928 to 0940 UTC (land to ocean 
return flight track). As expected, the values of opaque-channel brightness temperatures 
(Fig. 4a) increase with increasing offset from 183.3 GHz and are not responsive to the land- 
ocean background transitions. Relatively large brightness temperature changes are evident 
in the 89 GHz (—60 K) and 150 GHz (-18 K) window-channel observations (Fig. 4b) for both 
the island-ocean boundary at 9:09-9:13 UTC and at 9:28-9:32 UTC. The cooler brightness 
temperatures over ocean are the result of an increased contribution from the cold ocean- 
reflected cosmic background radiation. 
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A 355 nm Raman lidar [22] was deployed at Wallops Island during the five ER-2 
overflights and the ground-based experiments. The Raman component of molecular scatter- 
ing is characterized by a unique shift in frequency relative to that of the incident excitation, 
and thus unambiguously indicates the density of the scattering species. However, molecular 
Raman scattering is relatively weak compared to Rayleigh scattering from molecules or Mie 
scattering from atmospheric aerosols. Thus, the system could only be operated at night. At a 
given altitude, the water-vapor mixing ratio, defined as the mass of water vapor per unit 
mass of dry air [19], can be determined from the Raman-shifted return signals from water 
vapor and nitrogen. The lidar was deployed on a mobile platform, capable of making both 
zenith observations and range-height indicator (RHI) maps. A detailed description of the 
Raman lidar system and water-vapor retrieval algorithm can be found in Whiteman et. al 
[23]. A total of 40 Wallops lidar profiles were used in the clear-air comparisons. 

3.2 Ames Flights, May 1992 

The MIR flew on the ER-2 out of NASA Ames, Moffet Field, CA, during two flights on 
14-15 May and 15-16 May 1992 over the land and coastal waters of southwestern California 
(Table III). Statistics of MIR-observed brightness temperatures for these flights are shown 
in Table II. The two flight tracks were nearly identical and included overflights of San Nico- 
las and Santa Catalina Islands, and were within 60 km of Point Mugu Naval Air Station. 
The selection of MIR observations for cloud-free regions was facilitated by viewing data 
from the visible and IR channels of the Airborne Oceanographic Color Imager (AOCI Daed- 
alus), which was operating on both late -afternoon flights. For the clear-air comparisons, 
data from a total of six San Nicolas and Point Mugu radiosonde launches were used. The 
maximum temporal and spatial collocation differences between ER-2 positions and the 
radiosonde launches were 1 hour and 100 km, respectively. 

A nadiral brightness temperature time-series for the MIR is shown in Fig. 5 for the 
14-15 May flight from 0:12:00 to 0:30:00 UTC, corresponding to overflights of Los Angeles 
and its coastline (land to ocean) and Santa Catalina Island. The 183.3±1 GHz channel 
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(Fig. 5a) detects a gradual decrease in the water vapor of the upper atmospheric levels, 
manifested by a -15 K brightness temperature increase. Relatively large brightness temper- 
ature changes (40-80 K) are again evident in the window-channel observations (Fig. 5b) of 
both the land-ocean boundary at 0:19:30 and Santa Catalina Island at 0:22:30, particularly 
in the 89 and 150 GHz channels. Over land, brightness temperature fluctuations of up to 
20 K in the 89 GHz channel show the effects of the variable background emissivity. 

3.3 Radiosonde and Raman Lidar Atmospheric Profiles 

Radiosondes offer the most accepted means of obtaining verification measurements 
of the principle atmospheric profile parameters: temperature, pressure, and water-vapor 
density. For the Wallops experiments, single-, dual-, and triple-soundings were performed 
using AIR, VIZ, and Vaisala radiosondes. Balloons were launched at three-hour intervals 
from approximately 1800 to 0600 hours, local time. Multiple sounding packages were sus- 
pended from the same balloon, less than one meter apart, thus allowing each instrument to 
sample the same horizontal layer in the atmosphere. The three radiosonde packages differ 
primarily in the type of hygrometer used: the Vaisala type employs a capacitive hygrometer 
(Humicap) based on a hygroscopic dielectric material, and the AIR and VIZ types both 
employ carbon hygristors but use different relative humidity reporting (calibration) prac- 
tices. Data from a relative humidity triple-sounding is shown in Fig. 6, taken at 0700 UTC 
on 29 July at Wallops Island. Also in Fig. 6, a coincident Raman lidar profile is shown, 
observed from a site approximately 1 km from the balloon launch. Due to the weak Raman 
backscatter phenomena and sensitivity limitations in the lidar detectors, the maximum alti- 
tude for Raman -lidar derived water vapor profiles was limited to approximately 9.5 km. The 
radiosondes used in the Ames-based experiments were part of the Meteorological Rawin- 
sonde System based at San Nicolas and Point Mugu and employ Vaisala Humicap hygrome- 
ters. Temperature sensor differences among these radiosondes (AIR and Vaisala) of 1-2 K 
have been noted [10], and their absolute accuracy is ±0.5 K. A standard error-propagation 
analysis using IRTE calculations shows that the magnitude of the brightness temperature 


15 



bias associated with a ±1-2 K sensor bias at every level is 1 K or less for all MIR channels. 

The large relative humidity discrepancies evident in Fig. 6 between the radiosonde 
sensors is of concern in intercomparison studies [81. In particular, above 8 km the AIR sen- 
sor’s dry-end limit was found to be approximately 20% RH; this is consistent with previous 
hygristor studies [24]. The VIZ sensor, which uses a different calibration process than that 
for the AIR [25], and the Viasala sensor showed up to a 20% RH disagreement in the 
8-16 km levels. Here, the VIZ and Vaisala discrepancy is consistent with a previous inter- 
comparison of the two instruments [26]. Up to the 6 km level (500 mb) there is good agree- 
ment in AIR, Vaisala, and lidar water-vapor reporting. The relative humidity bias due to 
temperature sensor thermal-lag (-1-2K warmer than ambient up to the 500 mb level) [27] is 
reported to be -4% RH [25] for AIR and VIZ hygrometers. Apparently due to dry-end limit 
calibration difficulties the thermal-lag bias was not seen in the Wallops or Ames profiles and 
thus was not studied here. A standard error-propagation analysis using IRTE calculations 
shows that the brightness temperature variation associated with a hygrometer variation of 
+5% RH at every level is approximately +3.5 K and +2.0 K for the 89 and 150 GHz channels, 
respectively, and less than -1.5 K for the other channels. The implications of the upper-tro- 
pospheric relative humidity reporting errors on comparisons with MMW radiometric obser- 
vations are discussed in the next section. 

4.0 Clear Air Intercomparisons 

Radiative transfer calculations were performed using an iterative numerical model 
applicable for a horizontally-stratified clear atmosphere [28]. The Liebe water-vapor absorp- 
tion model [11] and the Rosenkranz oxygen absorption model [12] were used to relate the 
radiosonde and Raman lidar profile data to absorption coefficient profiles. Profile data were 
obtained by sampling the radiosonde measurements and lidar observations along a nearly 
uniform altitude grid from 0 to 20 km, i.e., -150 m spacing from 0 to 10 km, and 500 m spac- 
ing from 10 to 20 km. Of the four instrument types (AIR, VIZ, Vaisala, and lidar) the AIR 
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data were typically the most complete and included more than 60 levels. Some VIZ data suf- 
fered from grid spacings of up to 2 km below the 10 km level due to the data files becoming 
corrupted. In these cases, the VIZ profiles were augmented with AIR data for those levels. 
Also, a few VIZ and Vaisala data sets were empty above approximately 12 km and were aug- 
mented with standard mid-latitude exponential water-vapor atmospheres, 37° N for August 
at Wallops, and 33° N for May at Ames. For levels in the lidar data above 9.5 km, both AIR 
and Vaisala data were used to augment the profiles. For this study all radiosonde and lidar 
profiles were considered as zenith measurements. 

The effect of a non-ideal antenna gain pattern on the calculated brightness tempera- 
tures was significant for the zenith-directed ground-based observations. For example, using 
Eq. (16) with a representative antenna gain pattern, a brightness variation of +50 K from 
zenith to a scan-angle of ±45° (typical) biases the zenith observation by +0.7 K. For the air- 
craft-based observations, brightness variations were less than 10 K over the scan-range and 
therefore did not significantly perturb the nadiral measurements. 

4.1 Aircraft -based comparisons 

The results of the aircraft-based brightness temperature comparisons are shown in 
Tables IV and V for Wallops-based experiments over ocean and land backgrounds, respec- 
tively, and in Table VI for Ames-based experiments over an ocean background. The ocean 
was modeled as a specularly reflecting surface using the Fresnel reflectivities. No surface- 
foam was considered due to the relatively calm conditions throughout both experiment 
areas, and ocean surface roughness was neglected. Sea surface temperatures were deter- 
mined from seasonal means for the coastal Pacific in May and the Atlantic in July [29]. 
Land background emissivities were considered frequency independent and uniform with 
zenith angle 6 . 

The three comparisons with the radiosonde instruments in Table IV, and the one in 
both Tables V and VI, are generally good except for the 183.3±1 channel. Here the effects of 
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the AIR dry-atmosphere reporting errors are the most pronounced, and differences of 
-10.9 K are seen. Because of the sign (negative) of W„ (z, 6 ) at altitudes from 8-12 km 
(Fig. 2a), the excessive moisture reported in the higher levels by the AIR radiosonde reduces 
the computed opaque-channel (183.3±1 and 183.3±3 GHz) brightness temperatures. For 
example, a +20% RH perturbation from 8 to 20 km results in a brightness temperature per- 
turbation of -8 to -10 K at 183. 3±1 GHz. Similarly, VIZ comparisons also show a significant 
(but smaller) discrepancy for the 183.3±1 channel, which is also due to upper level hygristor 
reporting errors. However, since three of the nine VIZ profiles were augmented with AIR 
upper-level data, some of this discrepancy also can be ascribed to the AIR. Vaisala compari- 
sons exhibit the best agreement of the radiosondes types, with the 89 and 150 GHz channels 
showing the worst agreement for this instrument. For these two channels, the biases are 
nearly the same magnitude and are likely due to the under-reporting of water vapor in the 
lower atmosphere by the Vaisala’s Humicap sensor. 

The lidar comparisons exhibit their best agreement at 89 and 150 GHz (Tables IV 
and V). These channels are most sensitive to water-vapor profiles from 0 to 9.5 km, corre- 
sponding to the levels where lidar data is available. Lidar levels from 9.5-20 km were filled 
with the corresponding data from both AIR and Vaisala profiles. The lidar (w/AIR) upper- 
level comparisons, i.e., those for the remaining four channels, are similar to AIR results, 
particularly for the 183. 3±1 channel. Whereas this similarity is probably due to gross AIR 
dry-atmosphere reporting errors dominating the comparisons, the widely differing upper- 
level lidar(w/Vaisala) and Vaisala comparisons are more difficult to explain. Because the 
Vaisala hygrometer more accurately reports the upper-level water-vapor profiles, the discon- 
tinuous transition between the Vaisala and lidar profiles might be dominating the compari- 
son and causing the upper-level lidar (w/Vaisala) and Vaisala discrepancy. 

For the 89 and 150 GHz comparisons in Tables IV and VI, variations in ocean-sur- 
face emission can explain the standard deviations of up to -5 K. For example, a ±3% varia- 
tion in surface emissivity can lead to a 2 K and 1 K change in calculated 89 and 150 GHz 
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brightness temperatures, respectively. Also, the relatively large standard deviations m the 
lidar comparisons in Table V are likely due to modeling errors of the land-surface emissivity. 


4J2 Ground-based comparisons 

Brightness temperature comparisons for the two ground-based experiments at NASA 
Wallops are shown in Fig. 7 for 29-30 July, and Fig. 8 for 30-31 July. Local weather condi- 
tions at the outset of both experiments were similar and unusually clear and dry (-50-75% 
SRH) for coastal Virginia in July. During the first night’s experiment the atmosphere was 
characterized by an increase in relative humidity at lower levels, manifested by a bright- 
ness temperature increase (5-10 K) for the less opaque channels, 89, 150, and 220 GHz (Fig. 
7a). Conversely, for the 30 July experiment there was a decrease in relative humidity at 
lower levels, manifested as a 5-10 K decrease in the same channels (Fig. 8a). The sensitivity 
of these channels to changes in relative humidity is consistent with the behavior of the 

incremental weighting functions W p> (z) (Fig. 3). 

The 29-30 July brightness comparisons (Fig. 7b) were made over two time frames: a 
3 hour interval for AIR-based calculations, and a 1 hour interval for lidar- and Vaisala- 
based calculations. The AIR comparisons show a generally increasing discrepancy over the 
3 hour interval, indicating that the influence of that night’s increasing relative humidity is 
more pronounced in the calculations than in the MIR observations. The discrepancy change 
is largest for the 150 and 220 GHz comparisons, which is consistent with these channels’ 
higher sensitivity to relative humidity perturbations. The lidar and Vaisala comparisons are 
generally good over the 1 hour interval, except for the 150 GHz Vaisala comparison. This 
apparent bias may be due to a hygrometer induced error, given the 150 GHz channel is most 
sensitive to water-vapor perturbations. 

The 30-31 July comparisons (Fig. 8b) were done over a 5 hour interval for AIR-based 
calculations, 3 hour interval for the Vaisala-based calculations, and 1.5 hour interval for 
lidar-based calculations. For the AIR comparisons we see a smaller discrepancy change rela- 
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tive to the 29*30 July results, but still the largest discrepancy overall. Given the -5*10 K 
decrease in the observed window channel brightness temperatures over the evening, the 
AIR fails to indicate changes (drying) in the lower atmosphere. The Vaisala results show a 
negative-valued discrepancy, as was seen in the 29-30 July comparisons, but generally are 
in better agreement. As was found in the aircraft-based comparisons, this is most likely due 
to the Humicap’s under-reporting of water vapor in the lower atmosphere. The comparisons 
with the Vaisala do indicate the effects of the decreasing relative humidity over their 
3.0 hour interval. A decreasing discrepancy is seen in the lidar comparisons at 150, 220, 
183.317, and then 89 GHz. The 220 and 183.317 GHz lidar results are similar and the calcu- 
lations appear to over-report the effect of atmospheric drying relative to the MIR observa- 
tions. Since these channels have widely different sensitivities to water-vapor perturbations, 
their common behavior may indicate a possible MIR measurement bias for this time frame. 
However, the other two radiosonde comparisons, particularly the Vaisala results, do not 
indicate MIR measurement difficulties during this period. The 150 GHz channel appears to 
be most sensitive to the apparent bias in both the AIR and lidar calculations. 

The accuracy of these ground-based measurements is sensitive to the temperature of 
the MIR cold calibration target. The cold-target temperature is typically 250 K in-flight, and 
is near ambient temperature (-300 K) on the ground. The same noise or measurement 
uncertainty in the determination of the cold-target temperature, when the target is at 
300 K, produces a calibration error approximately four times greater than that when the 
target is at 250 K. That is, for a ±0.1% error in the cold-target temperature, the bias in a 
150 K calibrated brightness temperature would be approximately ±2.3 K for the 300 K cold 
target, and ±0.6 K for the 250 K target. This source of error alone, however, does not explain 
the disparate bias seen in comparisons based on observations made at the same, or similar, 
times, i.e., the AIR comparisons are almost all positive-valued, while the Vaisala are almost 
all negative-valued. 
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5.0 Conclusions 


The experiments described herein use the first atmospheric observations made with 
the wideband MIR channel set, and include unique observations at 220 GHz. The aircraft- 
based comparisons of computed-less-observed brightness temperatures at the two most 
opaque MIR channels (183±1 and ±3 GHz) show that the AIR dry-atmosphere relative 
humidity reporting errors produce discrepancies of -11 K and -6.5 K, respectively, over both 
land and ocean backgrounds. For the 89, 150, 183.3±7 and 220 GHz channels, and for the 
opaque channels compared against calculations based on other radiosonde types, the air- 
craft-based brightness temperature comparisons are in good agreement. Since different 
types of hygrometer-based radiosondes and reporting practices can introduce biases in 
upper-level water-vapor measurements, the potential effects of radiosonde inconsistency 
cannot be neglected in climatological and satellite calibration/validation studies. This fur- 
ther suggests that aircraft- or spacecraft-based passive microwave water-vapor profile 
retrievals can lead to more consistent measurements, albeit with reduced vertical resolu- 
tion. 


For the aircraft-based comparisons in this study, the most significant radiosonde- 
based sources of error are (in order of importance): 1) radiosonde hygrometer non-represen- 
tative calibrations and non-responsive behavior in dry environments, particularly for the 
AIR package; 2) hygrometer measurement uncertainty; 3) incomplete or corrupted water- 
vapor profile data; 4) temperature-induced humidity errors due to temperature sensor ther- 
mal-lag; and 5) inconsistent temperature sensor behavior among the representative pack- 
ages. Calculations based on Raman lidar humidity data are generally in excellent 
agreement with MIR observations of the lower troposphere. However, upper-level compari- 
sons are impeded by lack of lidar data above 9.5 km. The effects of atmospheric absorption- 
model uncer tain ties at millimeter wavelengths are not large enough to influence brightness 
temperature comparisons relative to the radiosonde errors. 

Ground-based window channel calculations (particularly at 89 and 150 GHz) are 
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more sensitive to hygrometer errors in the lower atmosphere than are airborne nadiral cal- 
culations. For example, the ground-based relative humidity incremental weighting func- 
tions for these channels are 5-20 times larger between 0-3 km than in the airborne case. 
Therefore, the effects of hygrometer uncertainties (±5 RH) on the calculations are consid- 
ered the dominant cause of the discrepancies reported here. Furthermore, differences 
between the hygrometer observations are evident by the disparate values for each respec- 
tive comparison, i.e., AIR-based calculations are positive-valued for each channel, and 
respective Vaisala-based calculations are negative-valued. Comparisons between MIR 
observations and lidar-based calculations show the best agreement with discrepancies typi- 
cally better than ±6 K. 
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Table I 

MIR Receiver Characteristics 


LO Frequency 
(GHz) 

IF (MHz) 

DSB Receiver 
Temperature (K) 

Nominal 
AT™ (K) 

89.0 

500-1500 

630 

0.13 

150.0 

500-1500 

860 

0.16 

183.31 

500-1500 

2000 

0.34 

183.31 

2000-4000 

2000 

0.28 

183.31 

6000-8000 

2000 

0.28 

220.0 

1000-4000 

2000 

0.26 
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Table n 

Summary of the observed nadiral MIR brightness temperatures used in the comparisons 

from Wallops and Ames flights 


Location (back- 

Frequency (GHz) 

ground type) 

89 

150 

183±1 183±3 

183±7 

220 


Wallops (ocean) 

(26 samples) 


Mean (K) 

220.1 

266.7 

Std. Dev. (K) 

10.41 

12.34 

Wallops (land) 

(8 samples) 


Mean (K) 

279.9 

286.4 

Std. Dev. (K) 

3.19 

2.51 

Ames (ocean) 

(8 samples) 


Mean (K) 

207.3 

252.4 

Std. Dev. (K) 

2.07 

3.27 


255.6 

270.3 

281.3 

282.1 

5.66 

5.13 

4.72 

5.03 

255.1 

270.0 

280.7 

285.7 

6.56 

5.42 

3.88 

2.61 

260.1 

275.3 

282.2 

277.0 

3.91 

2.68 

1.77 

2.48 


29 



Table m 

Summary of ER-2 Data Flights for Clear- Air Comparison Experiments (1992) 


ER-2 Flight # 

Date 

Time (UTC) 

Description 

92-089 

14-15 May 

2315-0440 

San Nicolas, Pt. Mugu 

92-090 

15-16 May 

2320-0445 

San Nicolas, Pt. Mugu 

92-131 

29 July 

0700-1325 

Flight track l a 

92-132 

30 July 

0707-1330 

Flight track l a 

92-134 

2 August 

0630-0920 

Flight track 2 b 

92-135 

3 August 

0300-0600 

Flight track 2 b 

92-140 

6 August 

0700-1310 

Flight track l a 


a. 2 passes over Wallops parallel to coast, 2 passes over Wallops perpendicular to 
coast, 1 round-trip to Cape Cod over ocean, 200 km off coastline. 

b. 2 passes over Wallops parallel to coast, 2 passes over Wallops perpendicular to 
coast. 
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Table IV 

Wallops Island aircraft-based comparisons over an ocean background of computed 
and observed brightness temperatures. Nadir values: computed - observed 


Instrument Type 


Frequency (GHz) 



89 150 

183±1 

183±3 

183±7 

220 

AIR Radiosonde 

Mean (K) 

(18 comparisons) 
3.42 1.43 

-10.90 

-6.47 

-2.98 

-0.432 

Std. Dev. (K) 

2.82 3.45 

3.61 

3.49 

2.27 

1.36 

VIZ Radiosonde 

Mean (K) 

(9 comparisons) 
3.62 0.50 

-6.33 

-2.19 

-0.41 

-0.462 

Std. Dev. (K) 

2.84 1.02 

1.59 

1.57 

0.94 

0.61 

Viasala Radiosonde 
Mean (K) 

(12 comparisons) 
-2.27 -2.49 

-2.09 

0.23 

1.04 

-0.033 

Std. Dev. (K) 

2.07 1.80 

3.06 

1.88 

1.34 

0.76 

Lidar (w/AIR) 

Mean (K) 

( 125 comparisons) 
-0.093 -0.41 

-11.10 

-4.60 

-1.54 

-0.282 

Std. Dev. (K) 

2.30 3.16 

3.20 

2.47 

1.89 

1.94 

Lidar (w/Vaisala) 

Mean (K) 

(69 comparisons) 
-0.141 -0.93 

-5.87 

-2.51 

-0.381 

-0.503 

Std. Dev. (K) 

1.80 1.83 

2.16 

1.43 

1.08 

0.86 


31 



Table V 

Wallops Island aircraft-based comparisons over a land background of computed and 
observed brightness temperatures. Nadir values: computed - observed 


Instrument Type 

Frequency (GHz) 

89 150 

183±1 

183±3 

183±7 

220 

AIR Radiosonde 
Mean (K) 
Std. Dev. (K) 

(5 comparisons) 
3.35 0.23 

0.89 0.54 

-10.40 

3.51 

-5.69 

2.75 

-2.34 

1.51 

-0.68 

1.04 

Lidar (w/AIR) 

Mean (K) 
Std. Dev. (K) 

(43 comparisons) 
-0.76 0.61 

19.24 4.77 

-10.90 

4.45 

-4.09 

3.10 

-0.67 

1.30 

1.73 

1.79 
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■Table VI 

Ames aircraft-based comparisons over an ocean background of computed and observed 
brightness temperatures. Nadir values: computed -observed 


Instrument Type 

Frequency (GHz) 

89 150 

183±1 

183±3 

183±7 

220 

Vaisala Radiosonde 
Mean (K) 
Std. Dev. (K) 

(16 comparisons) 
2.05 0.38 

4.02 5.63 

-3.35 

3.68 

-2.71 

2.37 

-1.18 

1.39 

-1.64 

3.13 
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7.0 Figures 

Figure 1: Gaseous absorption coefficient due to molecular absorption of oxygen and water- 
vapor, for (a) z=0 km, P=1000 mb, and p v = 17.0 kg/m 3 , (b) 2 km, 800 mb, 6.60 kg/m 3 , (c) 6 km, 
450 mb, 0.66 kg/m 3 , and (d) 12 km, 200 mb, 0.039 kg/m 3 . Specific frequencies of interest to 
this study are indicated by arrows. 

Figure 2: Airborne-based nadir-directed incremental relative humidity weighting func- 
tions for a 2-km scale-height exponential water-vapor profile over: (a) an ocean background, 
and (b) a land background. 

Figure 3: Ground-based zenith -directed incremental relative humidity weighting functions 
for a 2-km scale-height exponential water-vapor profile. 

Figure 4: MIR brightness temperature time-series for ER-2 flight 92-131 over the NASA 
Wallops Island facility: (a) 183.3±1, ±3, ±7 GHz channels, (b) 89, 150, and 220 GHz channels. 
The flight track is a west- and east-directed overflight of the island and Virginia coast at 
0900-0917 UTC and 0928-0940 UTC, respectively. 

Figure 5: MIR brightness temperature time-series for ER-2 flight 92-089: (a) 183. 3±1, ±3, 
±7 GHz channels, (b) 89, 150, 220 GHz channels. This flight track is an overflight of Los 
Angeles and San Nicolas Island for 0012-0030 UTC. 

Figure 6: Typical radiosonde relative humidity profiles from an AIR, VIZ, and Vaisala tri- 
ple-sounding and a Raman lidar profile. (Wallops Island Va., 29 July 1992, 0700 UTC.) 

Figure 7: Wallops ground-based experiment performed from 2200 to 0210 UTC on 

29- 30 July: (a) Time series of the six zenith-directed MIR observations, (b) Brightness tem- 
perature comparisons between MIR and AIR-, Vaisala-, and lidar-based calculations. 

Figure 8: Wallops ground-based experiment performed from 2200 to 0250 UTC on 

30- 31 July: (a) Time series of the six zenith-directed MIR observations, (b) Brightness tem- 
perature comparisons between MIR and AIR-, Vaisala-, and lidar-based calculations. 


34 





220 






Altitude (km) 



W*fe)(Kknr'%') 

Figure 3 


37 












Altitude (km) 



Relative Humidity (%) 

Figure 6 


40 


-Observed Brightness (K) Observed Brightness (K) 






-Observed Brightness (K) Observed Brightness (K) 



22:00 23:00 24:00 1:00 2:00 3:00 

Time (h:mm, UTC) 



22:00 23:00 24:00 1:00 2:00 3:00 

Time (h:mm, UTC) 

Figure 8 


42 




OPTIMAL CALIBRATION OF RADIOMETERS USING WIENER FILTERING 
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Abstract 

The application of Wiener filters to the problem of ra- 
diometer calibration is investigated. The calibration pro- 
cess uses voltage measurements of two targets of known 
temperature and emissi vity, to derive a gain and offset for 
w h scan of the instrument. The actual measurements, 
however, contain integration noise which limits the accu- 
racy of the calibration. This research focuses upon the 
application of statistical Wiener filtering theory to de- 
velop optimal (minimum mean squared error) linear es- 
timates for the gain and offset. The filters are tested for 
data from the Millimeter- Wave Imaging Radiometer. A 
significant reduction in calibration noise relative to that 
of a conventional infinite impulse response filter is demon- 
strated. Since the filter was derived for the conditions of 
statistical stationarity, a method for detecting jumps in 
the gains and offsets was developed. 

1. Introduction 

In microwave radiometry, measurements of hot and 
cold blackbody targets of known temperature and emis- 
sivity typically are used for calibration. The technique 
involves estimation of two slowly time varying processes 
(the gain 01 and the offset b) from noisy observations of 
these targets. Typically, these targets are observed once 
per scan producing voltage output v — mT^+b+n where 
a represents additive white Gaussian noise and Ta is the 
antenna temperature. Using observations of the hot and 
cold targets (v#, and vci) coupled with the known an- 
tenna temperature, the single scan estimates for the cal- 
ibration quantities are: 

VBi ~ VCi r _ v CiTm — VHiTci f jx 

*"* ~ Tui - Ta 1 ~ Tb* ~ Ta 

However, because of the presence of the noise a, these 
estimates are necessarily inaccurate. 

Typically, this error is reduced by applying a sin- 
gle pole infinite impulse response (1IR) filter. This filter, 
however, is slow to adjust to shifts in the gain and off- 
set statistics experienced during normal operation of the 


radiometer. A better filter would make use of the statis- 
tics of the gun and offset shifts to optimally reduce the 
noise component of the signal by incorporating many sin- 
gle scan observations into the development of the desired 
single spot quantity. 

We begin by defining: 

rhi = ^2 fhi-k = X* “ft **'“* ( 3 ) 

ft * 

where the and u\ are the filter weights for the respec- 
tive process m and b. When these weights are chosen to 
minimize the mean squared error of the the estimate, this 
b the discrete Wiener filter. Thu paper explores the ap- 
plication of Wiener filtering to the problem of radiometer 
calibration. 

2. Discrete Wiener Filter Theory 

The Wiener filter coefficients can be obtained by 
fninimi»tng the mean squared error between the filtered 
estimates and the actual desired quantities. Letting zj 
denote the desired variables (i.e. true gun or offset) and 
t n denote its single scan estimate the error b defined as: 

ti = *i —*i 

N 

= (z<- 53 “ft^-ft) ( 3 ) 

km-N 

where a 2 N + 1-point function b used. Since process- 
ing b performed off-line, the filter can be non- causal. A 
non-causa) filter b used in recognition of the fact that 
future statbtical trends are relevant to the calculation of 
the current quantities. If a causal filter were used, all in- 
formation about the future signal would be lost and the 
error of the system would increase. 

The mean squued error b defined as: 

c\ = Eiixi-ztf] 

N 

= £[(*<- 5Z (*) 

ft s-N 



Using the autocorrelation function estimates m and 
h, the filter 1 5 was computed: 

IF sr (Ret + *Rnn) *^rv 0^) 

To Insure unbiased estimates of m and b, unitarity (i.e. 
£ic= 1) was enforced by scaling the filter coefficients. 



The physical interpretation of the impulse at t = 0 
in TL* is the variance of the high frequency fluctuation 
of the gun and offset processes due to integration noise. 
In the absense of gain and offset fluctuation noise, the 
n'inim nm standard deviation of a constant temperature 
scene is [2, pp. 358 - 368]: 

A7*ms = 


where B is the bandwidth of the the channel, r is the in- 
tegration time for the instrument, and T,y $ is the system 
temperature. This equation defines the relative sensitiv- 
ity of the instrument. As gain and offset noise increases, 
additional terms must be included to account for the fluc- 
tuations increasing the minimum sensitivity to 

AThms = T.„ + + (13) 

^ grading the performance of the radiometer. 

The improvement of the A Trms using the Wiener 
filtering technique was demostrated by processing dear- 
air radiometric data over ocean for which the scene 
brightness temperatures are nearly constant. Prior to ap- 
plying the filter to data, the estimated parameters from 
the autocorrelation functions were used for theoretical 
comparisons to test the accuracy of this model. From 
Eq. 1, the single-scan gain fluctuations are: 


nu 


t ’Hi ~ t >Ci 

Tm — Ta 


Th i - Ta 

__ VHA - VCA nH_-nc 
~ Tm — Ta Tm - Ta 


= 171 Ai + 


njf -nc 

Th, — Ta 


where m Ai is the true quantity (i.e. no noise). Substi- 
tuting into Eq. 4, a 7 , can be predicted by: 


_2 _ 


n h - »»c 


= *tt"u -("* + ? 

__ 2((m)ATRMs) 7 
~ (Tm - Ta? 


This allows comparison of the measured height of the 
impulse function to a theoretical calculation. 


For these comparisons, we use a modified version 
of Eq. 12 which accounts for the RC integration and 
oversampling technique [7] used in the MIR. 


ATrms 

P(N L ,p) 



where t = N L At, P(N l ,p) is a derating factor with 
Si total samples during calibration stares, and p — 
cxp(— &r) is the correlation between subsequent calibra- 
tion measurements. The actual sampling interval is At. 

As At decreases relative to the RC time constant, con- 
secutive voltage measurements contain information from 
prior samples still held by the integration filter. Since 
the previous signals have not yet decayed, new measure- 
ments contain less additional information than the com- 
pletely uncorrelated case (e.g. Ai » RC). The P-factor, 
therefore, corrects the ATrms in Eq. 12 due to this cor- 
relation. For the MIR, P = 1.48 with N L = 27, RC = 
6.8 ms and a sampling interval of approximately 6.1 ms. 

For the following noise calculations (Table 1), Tm — 
Ta — 100 K. By definition, r = N,T, (2] where N, is 
the number of samples taken at a specific location (spot) 
in the rotation of the mirror (AT, = 13) and T, is the 
sampling interval (6.1ms) yielding a value r » 80ms. It 
is noted that the three window channels (89, 150, 220 
GHi) all have good agreement with the theoretical cal- 
culation and the estimate obtained from the autocorre- 
lation function. The three 183 GHr channels, however, 
have greater values for the estimated noise than the the- 
oretical prediction (5 times greater). This does not imply 
failure of this method, but rather implies the presence of 
a noise source greater than the theoretical limit (i.e. Eq. 
13 rather than Eq. 12). Calculations from other flights 
demonstrate similar numerical results (i.e. good agree- 
ment for window channels end similar difference factors 
for the 163 GHt channels). 

The Wiener calibration filter is applied to the gain 
and offset sequences and the brightness temperatures are 
subsequently computed. The Wiener-filtered data is then 
compared with the IIR-filtered data for evaluation pur- 
poses using data collected over clear air regions of flight. 
Clear air regions were located using coincident measure- 
ments with the lidar. Since the lidar measurements verify 
clear air at nadir positions, the MIR comparisons have 
been limited to these spots to eliminate the possibility of 
atmospheric structure creating increases in the variance 
of the scene temperature. However, full scan calculations 
have been performed with results that are consistent with 
the nadir spots in areas believed to contain clear air. 

The standard deviation of the dear air measure- 
ments is calculated using the along track nadir spots for 
both the I1R and Wiener filtered case. In order to obtain 
statistically significant comparisons, only regions with a 


Table 1: Comparison of Estimated and Calculated Noise 
Variance 


Channel 

ATrM5 

<n 

— 

89 

0.199 

2.24e-7 

2.5e-7 

150 

0.245 

3.35e-7 

4e-7 

183 ± 1 

0.420 

1.58e-6 

8.15e-6 

183 ±3 

0.297 

8.15e-7 

4.25e-6 

183 ±7 

0.297 

8.31e-7 

4.5e-6 

220 

0.323 

2.4e-6 

3.4e-6 


Table 2: HR versus Wiener Filter A Trms Comparisons 
Scene 1 


Channel 

A7)tws 

IIR 

ATrju5 

Wiener 

89 

1.1195 

1.0738 

150 

0.4221 

0.3544 

183 ± 1 

1.7091 

1.3039 

183 ± 3 

1.1591 

0.9134 

183 ± 7 

1.6537 

0.6108 

220 

0.8920 

0.5239 


Table 3: HR versus Wiener Filter A Trms Comparisons 
Scene 2 


Channel 

ATrats 

IIR 

A7/LMS 

Wiener 

89 

0.9973 

0.9837 

150 

0.4307 

0.3187 

183 ± 1 

2.2003 

1.2264 

183 ± 3 

1.5338 

0.7959 

183 ± 7 

1.1878 

0.8391 

220 

0.9499 

0.5446 


Table 4: Jump Detector Performance 


Jumps Locations 

Jump Height 

Jump Detection 

700 

0.7078 


1200 

-0.7078 


1500 

-1.4138 

X 


1.4138 

X 

2100 

1.4734 

X 


-1.4734 

X 

2700 

0.4262 



-0.4262 


3300 

1.2315 

X 

3800 

-1.2315 

X 










Figure 5: Simulated eequeoce for testing jump detector. 



Figure 6: Simulated sequence with jumps located as de- 
picted in Table 1. 



Figure 7: Z-ratio test results. Notice the impulse func- 
tions at the jump locations. 
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Abstract 

Partially polarimetric measurements of thermal emission from a striated 
water surface at 91.65 GHz illustrate the potential for remote sensing of water wave 
direction by passive microwave radiometry. The three Stokes parameter measure- 
ments were made using a precision polarimetric radiometer trained on a rotatable 
water wave tank at several elevation angles from near nadir to near gracing. The 
polarimetric measurements are well corroborated by calculations using a tilted-facet 
geometrical optics model for the water surface emission and scattering. Multiple 
scattering of the incident background radiation is incorporated for observation an- 
gles approaching grazing. The downwelling background brightness is computed using 
an atmospheric radiative transfer model. We show that azimuthal brightness varia- 
tions in the third Stokes parameter are in phase quadrature with the first and second 
modified Stokes parameters. For observation angles near ~ 60 - 70° from nadir the 
first three parameters have particularly large azimuthal brightness variations, and 
thus have significant potential for measuring ocean wave direction. Moreover, the 
azimuthal brightness variations caused by water waves are not negligible for many 
passive microwave atmospheric sounding and surface remote sensing purposes, even 
at nadir. A range of elevation angles resulting in minimal azimuthal variations is 
identified. 
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1. Introduction 

The utility of vertically and horizontally polarized brightness temperatures 
in passive microwave remote sensing of terrestrial atmospheric and surface proper- 
ties is recognized by many investigators. Studies of dual-polarized thermal radiation 
have benefited from satellite data using instruments such as the Scanning Multi- 
channel Microwave Radiometer (SMMR) [Njoku et al ., 1980] and the Special Sensor 
Microwave/Imager (SSM/I) on the DMSP Block-5 platforms [Hollinger et al., 1990]. 
Here, the large contrast between T v and 7* (the first two modified Stokes 1 parameters, 
as defined in the Earth’s “natural” polarization basis) has been used to improve mea- 
surements of ocean surface wind speed, atmospheric clouds, rain and water vapor, and 
facilitate sea ice detection. Recent theoretical and experimental investigations have 
also shown measurable differences between T v and 7* caused by oriented atmospheric 
ice particles [Evans and Vivekanandan, 1990; Adler et al., 1990]. 

Many techniques for the measurement and interpretation of the first two pa- 
rameters have been published [e.g., Stogryn, 1967; Spencer, 1986]. However, compara- 
tively few experimental studies of the utility of the third and fourth Stokes parameters 
( Tu and Ty) for Earth remote sensing have been made. Primary among these are 
airborne linearly-polarized brightness measurements over open wave-covered water by 
Dzura et al. [1992]. These measurements showed signatures in T v , 7/,, and (indirectly) 
Tu that were related to the wind vector. In other experiments, values of Tu as high as 
~ 50 K have been measured over rough periodic moist soil at a frequency of 10 GHz 
[Veysoglu, 1991; Veysoglu et al., 1991; Nghiem et al., 1991]. The measurements were 
corroborated by an electromagnetic diffraction model based on the method of mo- 
ments. Currently, there is increasing interest in passive microwave remote sensing of 
land and ocean parameters using fully-polarimetric measurements. 

In this paper, we present the results of experiments investigating the rela- 
tionship between the first three Stokes 1 parameters of the upwelling radiation field 
over a striated water surface and the associated water wave parameters. Direct mea- 
surements of T v , 7\ and Tu were made over a rotatable outdoor water wave tank 
using a well-calibrated 91.65-GHz polarization correlating radiometer developed by 
the authors [Gasiewski, 1990; Gasiewski and Kunkee, 1992, 1993]. A polarized cali- 
bration load was used to achieve an absolute calibration accuracy of approximately 
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2 K for each Stokes parameter. Standing gravity waves were generated in the tank; 
by rotating the tank the water wave vectors could be oriented at arbitrary azimuthal 
angles with respect to the radiometer’s plane of observation. 

The polarimetric measurements were corroborated by a geometrical optics 
(GO) model of the water surface scattering and emission. The GO model considers 
the water surface to be composed of specularly-reflecting facets, with the upwelling 
radiation field computed as a sum of contributions from each facet [e.g., Ilyin and 
Raizer , 1992]. The model accounts for multiple scattering along ray paths, and yields 
values of the full Stokes vector at observation angles from normal to near grazing. 
Estimates of the downwelling background radiation field are obtained from radiosonde 
data and using the numerical radiative transfer model of Gasiewski and Staelin [1990]. 

Computed GO and measured results are presented to illustrate the rela- 
tionship between the first three upwelling Stokes’ parameters and the water wave 
direction. The peak-to-peak amplitude and phase variations of T v , Th and Ty with 
azimuthal angle are shown to depend upon the water wave slope distribution, scan 
elevation angle and the atmospheric state. A numerical analysis illustrating the sen- 
sitivity of the upwelling radiation to various atmospheric and surface parameters is 
presented. Finally, implications of the investigation with regard to passive microwave 
remote sensing of ocean wave direction and passive atmospheric sounding are dis- 
cussed. 

2. Experiment Description 

The polarization-correlating radiometer was implemented by modifying a 
91.65-GHz duai-linearly polarized radiometer (Fig. 1). The va and vb channels are the 
analog outputs of a conventional dual-polarized Dicke radiometer. The vc channel is a 
simple adding correlator with a post-detection summing circuit to cancel the relatively 
large orthogonal-mode signals [Gastews&t and Kvnkee, 1993]. The analog signals cor- 
respond to each of three radiometric parameters: the two orthogonal feedhorn mode 
powers v A oc T v and vb oc Th, and a third parameter vc oc cos 7 (A<p)Ty + sin 2 (Av?)Tv. 
Here, A ip is a variable phase shift added to one of the two local oscillator signals. 
The vc channel is a linear combination of the in-phase and quadrature correlations 
between the orthogonal-mode field amplitudes. By adjusting A<p, the third channel 
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can be made sensitive to either Ty or Tv- The total electrical path lengths of the 
two orthogonal mode channels were equalized to well within one correlation length 
l c = v v {W ~ 5 cm by adjustment of A/, where v p ~ 10® m/sec is the IF phase velocity 
and W = 2 GHz is the IF bandwidth. 


The water wave observations were made over a square rotatable wave tank 
measuring 1.2 m on a side and filled to a mean depth of ~ 7 cm (Fig. 2). The 
tank size was large enough to subtend the radiometer beam to approximately the 
third pattern null, thus minimizing beam spillover. A 1.2-m long vertically-oscillating 
plunger coupled by eccentrics to a variable speed motor excited standing water waves 
with peak amplitudes from 0.3 to 0.5 cm. The dispersion relation foi small-amplitude 
water waves is [ Phillips , 1980] 

u 2 = gk 

where p is the water density, t is the surface tension, d is the water depth, w is the 
angular frequency and k is the surface wavenumber. Accordingly, the plunger was 
driven between 1.4 - 1.9 Hz to provide wavelengths adjustable from ~ 7 to 12 cm. The 
wave profiles were approximately sinusoidal, and the 3-dB footprint of the radiometer 
subtended several water wavelengths. 

All experiments were performed outdoors under clear-sky background condi- 
tions. The radiometer’s antenna beam was trained on the wave-covered surface using 
flat subreflectors. Care was taken to insure that a cold sky background illuminated 
the water surface from all specular facet angles. Due to some minor feedhorn align- 
ment problems the polarization basis of the instrument was skewed 15° relative to 
the v — h basis of the water surface. This alignment error was removed during data 
processing using a Stokes’ parameter rotation as described in Gasiewski and Kunkee 

[1993]. 

To minimize measurement errors caused by slightly imperfect sinusoidal 
water waves (e.g., waves exhibiting small travelling-wave components, reflected sec- 
ondary wave components and small amplitude inhomogeneities) at least one full rota- 
tion of the tank was performed. Coincident radiosonde measurements of atmospheric 
temperature and water vapor profiles were obtained for computing sky background 
brightness temperatures. All radiosondes were launched from Athens, GA, located 


i + £ 

P9 


tanh(fcd) 


( 1 ) 



100 km E-NE of the experiment location. 


3. Geometrical Optics Model 


Within the GO model, a non-flat surface is modelled by a distribution of 
specularly-reflecting facets, each of which contributes to the overall upwelling radi- 
ation in accordance with the Fresnel reflectivity relations [Tsan^, et a/., 1985]. The 
surface material is assumed to be homogeneous and isotropic with dielectric param- 
eters and illuminated by a clear-air background radiation field from a medium 

with parameters = e 0 , pj = fi c . 


Consider plane wave reflections into the observation direction k, from the 
k th facet (Fig. 3). The incident wave direction is: 

h k = (7- 2n*n fc ) k t (2) 

where the fi* is the facet normal vector and / is the identity matrix. Let the illumi- 
nating wave from this direction be 

Ei k = ( E Vik V{ k -+■ Eh, k hi k ) e J ** (3) 


where the incident vertical and horizontal basis vectors and hi k are defined ac- 


cording to: 


v tk = h tk x k tk 


? k, k x z 
h ik = — 


\k, k x z | 

The contribution to the scattered field from this facet is: 


(4) 


E. k = (E v , k v t + E h , k h,)e jk ‘ 7 


(5) 


where the scattered field components are related to the illuminating wave components 
by the scattering function matrix: 





( 6 ) 


and 

r t t kg x z ... 

v t = h t x kg hg — -x — • (7) 

|*j x z | 

The scattering function matrix elements are computed using the Fresnel reflectivity 
relations for the TE and TM wave components in the frame of the facet: 
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where 


/w* = (Si* ' Pi k ) Rv{ 9 i k ) (S. • Pn ) + ($» ’ 9 k) Rh(Oi k ) (*>* • ?*) 

fvh k = (h ik * Pi k ) Rv(Oi k ) (t>. * p,» ) + (h ik * 9 k) Rh(Oi k ) ( 0 , • q k ) 

hv k = ( v ik • Pi k )Rv{0l k ) (*». * Pn) + (*>•* * 9*) R h( 6 ‘ k ) (*• * 9k) 

fhh k = ( h, k ■ Pi k ) Rv( 0 l k ) (h, • P n ) + (*i* • 9 k) Rh{ 0 i k ) (h t • q k ) 

and the Fresnel reflectivities are: 

i?i cos 6 lk - Tji Jl - (&) sin 2 6 lk 

R,(e lk ) = / \ k — 

Tft cos &i k + i? 2 yl ~ (£) sin 2 0/ 4 

Tfrcos 9 tk - 9i^l ~ (^) 2 sin 2 ^ 


Rh(6i k ) — 


Vicos 9 lk + 171 - (fe) 2sin2 ^* 


( 9 ) 

(10) 

( 11 ) 

( 12 ) 


( 13 ) 


( 14 ) 


where ifa and 772 2-re the intrinsic impedances and k k and &2 are wave numbers 
of free space and the surface material, respectively. Here, $i k = cos _1 (n* • k,) is the 
angle between the facet and the observation direction. The facet unit vectors are: 

_ k, x n k _ kj x n k 

qk |fc, x n*| \ki x n*| 

Pi k = q k x h k (I 5 ) 

p» k — q k x k t 


The scattered radiation T tk (ic t ) from a facet can be related to the incident 
radiation from direction k, k using the Stokes’ matrix: 


Ti k (k,) 


T., k 

T tkk 

9»u» 


1/v.J 2 

\u k \ 2 

R*{fvv k fvh k } 

-hn{fvvkflh k } 


IA..I* 

1/i.h* 1 2 

R*{fhv h fhhk ) 

—lm{fh%> k fhh k } 


2 Re{fw k r hvk ) 

2Re{fvk k f; Hk ) 

+ /vh k /ht, k } 

“I m {fvv k fhh k “ fvh k fhv k ) 


. 2I*n{fw k fh Vk } 

2Im{/vk»/hh»} 

+ Afc k /hi> k } 

R e {fvv k fhh k ~~ fhv k } . 

L 


where T IJt is the incident Stokes’ vector for direction k ik . For a single-scattering facet 
Ti k is the unpolarized downwelling brightness field (i.e., Ti yk = Ti hJt and Ti Uk = T, v * = 


( 16 ) 


_H JS 
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0). All facets are assumed to scatter incoherently. Thus, their contributions to the 
total scattered radiation field T(k,) add: 

T.(k.) = r r T. t (k.)p(aJ)Jad0 (17) 

In the expression above, p( a, (3) is the probability of a facet having slopes a = df /d x 
and /? = df /dy where the surface is described by z = f(x,y). Accordingly, 

^ —ox - fiy + z 

” y/Sr+WT l 



For a single-scattering facet, the thermal emission T ek (k g ) can be polarized, 
but Tv is zero. To see this, consider a facet of kinetic temperature T. Using detailed 
balancing, the principles of reciprocity and conservation of energy can be used to find 
the thermal emission. In the p, k — q tk basis this is: 



(l-liWjnri 

o 

o 


(19) 


Using the rotational transformation, this is written in the v — h basis as: 




(1 - cos 1 *,!*.(», ,)| ! - rin J *,!*»(«, ,)| J )r 
(1 - sin 1 Wfl v (0,,)| 1 - cos 1 Wfl»(«,J| 1 )r 

sin2#.(|/l.(«,.)| 1 -|i! t («,.)| 1 )r 

0 


( 20 ) 


where <p T is the angle of rotation between the two bases: 

<t>T = 77 — f 7 cos (kg • q k ) 

|n* • 


( 21 ) 


Thus, facet emission with nonzero Tjj is possible, but 7V emission is always zero 
provided that the medium is isotropic. As in the case of the scattered radiation the 
orientationally-averaged emission becomes: 


T'(k.) = r r T ek (kg)p(a,P)dadp 

J — OO •/— OO 

The total upwelling radiation field is: 

T(k.) = T.(%.) + T.fy 


( 22 ) 


(23) 
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For deterministic surfaces, p(a,/3) is replaced by the fractional projected 
area of a given surface element: 


p{a,p)dad@ 


(n t 

]J^h)dS 


(24) 


Here, (n* • k , ) accounts for the projected area of the facet in the scattering direction, 
and the denominator normalizes the function. The surface integral is over only the 
portion of the surface that can be seen from direction % t , thereby accounting for 
shadowing. Written as an integral over the x — y plane, the upwelling radiation 
becomes: 


T(h) 


f°° f°° PVfc&) + T $k (k.)] ( n* • k,)J 1 + a 2 + Pdxdy 

J-oc J — QO ’ 

/ / (rc* • k t ) Jl + a 2 + P 2 dxdy 

J-oc J-oo 


(25) 


where the Jacobian of the transformation from surface coordinates to the x — y plane 
is used. Again, shadowed regions are excluded. 


Consider the GO model applied to the sinusoidal surface 
f{ x iV) = Asin(2;rx/A), where A is the peak wave amplitude and A is the surface 
wavelength. Here, a = (27rA/A)cos(27rx/A) and fi = 0. At large enough observation 
angles, multiple scattering either with or without shadowing can occur. Thresholds 
for the onsets of both multiple scattering and shadowing, in terms of 6 t , are shown for 
sinusoidal surfaces as a function of wave height for <j>, = 0 in Fig. 4. For 6, beyond 
the multiple scattering threshold, the k, k associated with some surface facets intersect 
the surface at other facets (Fig. 5). Accordingly, the unpolarized downwelling back- 
ground brightness T ik , must be replaced by the brightness of the facet at the point 
of intersection: 

r.. = t£> (£,„) + t*; 1 (i.j (26) 

where the superscript < 2) indicates radiation from the second facet. The process 
is repeated if the k, k for the second facet intersects the surface again. Note that 
the incident brightness illuminating the primary facet is not necessarily unpolarized. 
Consequently, Tv is not necessarily zero in the case of multiple scattering, even if 
the surface medium is homogeneous and isotropic. Shadowing occurs when the k t 
from any facet intersects the surface; shadowing is always accompanied by multiple 
scattering. For example, for A/A = 0.05 multiple scattering can occur for 6, ~ 51° 
and shadowing can occur for for 6, > tan -1 [A/(27rA)] rs 72°. 
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Use of the GO model is justified by considering the domain of applicability 
of the Kirchoff approximation (KA). As discussed by Wirgin [1983], the KA for a 
sinusoidal surface is applicable for hA/A 2 < 0.011 for 6, = 0, where A is the electrical 
wavelength. Applying this criterion to a periodic water surface with A ~ 10 cm and 
h ~ 0.5 cm, it is seen that the KA is indeed applicable at 91.65 GHz (A = 3.3 mm) 
for normal observation. As shown in the next section by the good agreement be- 
tween theory and experiment, this KA criterion appears conservative for polarimetric 
emission studies. For surfaces that are uniformly illuminated over at least several 
water periods A, (20) need be performed over only a single period. That is, any local 
brightness variations over the surface are adequately filtered by the antenna beam. 
This is the case for the experiments described herein. 

4. Computed and Measured Results 

Calculations of T v , T h and Ty based on the GO model for a sinusoidally- 
striated fresh water surface are shown in Fig. 6a-c. The calculations use h/A = 0.05 
and axe for observation angles 0 < 6, < 68° and 0 < 4>, < 90°. The brightnesses for all 
other quadrants in <f> , can be obtained by symmetry: T v and Th have even symmetry 
about <f>, = 0° and (f>, = 90° [i.e., T v {<f> t ) — T v (—<f> t ) and T v (90 -I- (f > t ) = T v (90 — <£»)], 
while Tjj has odd symmetry about <j> t = 0° and <f> , = 90° [i.e., Ty(<f> g ) = —Ty(—<f> 3 ) 
and 71/(90 + <f> t ) = -7V(90 - <f> 3 )}. The dielectric parameterization of Klein and 
Swift [1977] at 91.65 GHz and at a temperature of T w = 290 K is used for the 
water permittivity, resulting in £2 = 7.39 — j 12.38. The unpolarized downwelling 
radiation field is computed using the radiative transfer model described by Gasiewski 
and Staelin [1990] with a US standard atmosphere temperature profile interpolated 
to July at 34° N latitude. The water vapor profile was exponentially decaying with 
a 2 km scale height and 50% surface relative humidity (SRH). The figures for T v and 
Th show only the differences with respect to the vertical and horizontal brightnesses 
of a calm water surface, indicated by T v (e) and T[ c) (respectively). 

Three characteristic ranges in the relationships between T v , T h and Ty and 
water wave angle can be identified, each corresponding to particular ranges of 6,. 
When 0° < 6, ~ 35°, T v decreases, 7* increases and Ty remains negative as 4>, 
is moved from 0° (parallel to the water wave vector) to 90° (perpendicular to this 
vector). At nadir, (6, = 0°) purely sinusoidal variations of ~ 2 — 4 K are found in all 
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three parameters; this follows from elementary properties of the Stokes’ parameter 
rotation matrix. At 0 t ~ 40°, T v , and Tv are essentially constant in 4> t \ such angles 
would not be useful for observing wave direction. When 0 , ~ 45°, T v increases and 
T v remains positive as 4> t is moved from 0 to 90°. However, for 0 , ~ 50°, multiple 
scattering causes non-monotonic behavior in 7* as <f> M varies. 

The above behavior can be conveniently described using the azimuthal Fourier 
amplitudes B an and phases $ on , defined by: 

T a {0„ 4 > ,) = 53 “ T o C) ( 27 ) 

n=0 

where is the calm water brightness. In general these coefficients are functions of 
0„ as illustrated in Figs. 7a-c for n = 0,2, The B a0 (Fig. 7a) are the azimuthally- 
averaged brightness changes caused by the water waves. Even for this DC case 
values of B v o and Bh o exceeding 5 K at nadir and 20 K for T/, at 0 t ~ 55° are seen. 
Thus, the presence of small water waves can significantly change the azimuthally- 
averaged brightness. For 0 t ~ 50° the change is positive; this is a result of an 
increased fraction of the total scattered radiation originating from lower (and hence 
warmer) sky elevation angles. However, for 0 t ~ 50° this change is negative for 
T v and trending toward negative for 7*. For T v this is a result of a reduction in 
surface emission caused by the spread in facet angle distribution, the well-known 
null in reflectivity near the Brewster angle is not as distinct nor deep when waves 
are present. For the decreasing trend is a result of multiple scattering, more of 
the radiation is multiply-scattered and originates from higher (and hence colder) sky 
elevation angles when waves are present. As expected from symmetry, Buo = 0- 

The B q2 (Figs. 7b and c) and associated phases clearly show the three 
characteristic ranges of 0 t . Over 0° 5: 0* ~ 35°, B v 2 and B ^2 are of comparable 
magnitude but opposite sign. At nadir ( 0 , = 0°), these two coefficients are precisely 
equal and opposite (B v 2 — Bm, $v2 = $/»2 = — 180°) and = 2 * |-H„ 2 | — 

2 * \B h2 \] at nadir all B on are identically zero except for n = 0,2. Near 0, ~ 40° 
all B a 2 exhibit nulls and sign reversals. A second sign reversal is observed in Bh7 at 
~ 58° as multiple scattering becomes significant. Note that the $ 1/2 is always ±90°, 
while $> u2 and are always either 0° or 180°. That is, T v is in phase quadrature 
with T v and 7\. One important consequence of this is that observations of Tv along 
with coincident observations of T v and T h will facilitate the retrieval of wave direction 
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by removing ambiguities in the quadrant of the wave azimuthal angle [ Kunkee and 
Gasiewski, 1993]. All B an are zero for odd n, as expected for a symmetric surface. 
Higher order coefficients ( B an for n > 2) are generally nonzero, but usually too small 
to be of practical importance. An exception is B which has a maximum amplitude 
of ~ 4 K at 6, fa 58°. 

Measurements of T v , 7*, and Tv at B, = 20°, 39° and 65° using the 91.65 
GHz polarimetric radiometer are compared with computed GO brightnesses in Figs. 8- 
10. These three observation angles lie within each of the characteristic ranges de- 
scribed above. All measurements were over a sinusoidally striated water surface with 
h ~ 0.5 cm and A ss 10 cm under clear skies. The radiometric noise on each measure- 
ment is ~ 0.2 K and the bias is ~ 3 — 5 K. To minimize the effects of instrument drift, 
the data have been corrected by first subtracting the measured calm water brightness 
at each azimuthal angle, then adding the average measured calm water bright- 
ness (T^). The calm water measurement was made after turning off the wave tank 
plunger and allowing the waves to dissipate. 

Overall, the amplitudes and phases of the azimuthal variations in the mea- 
sured brightnesses corroborate the GO model well. As predicted, the amplitude vari- 
ations at 6 t = 39° are relatively small, while near nadiral ( B t = 20°) and near-grazing 
( 6 , = 65°) measurements both show significant variations. In general, these varia- 
tions are large enough to warrant consideration in many remote sensing problems. 
The sharp features in 7), and Tv for 6, = 65° are the result of multiple scattering 
of order two (double scattering). Indeed, due to the larger local reflectivity for h— 
compared to v — polarization multiple scattering is more strongly manifested in 7* 
than in T v . The dashed lines in Figs. 10a and b show GO simulations obtained when 
multiple scattering is neglected. In this case the model seriously conflicts with the 
measurements, thus clearly demonstrating the need to account for multiple scatter- 
ing at large observation angles. The absence of sharp features in the measured 7\ 
and Tv data at 6, = 65° is primarily the result of imperfect sinusoidal water waves: 
small capillary waves from air gusts and vibrations were often superimposed on the 
dominant gravity waves. 

Other discrepancies between the measured and computed brightnesses in- 
clude small biases of ~ 3 K in the measured values of T v for B t = 20° and ~ — 4 K in 
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Tf i for 6, = 39°. In addition, the amplitude of Ty for 6, = 20° is only ~ 60% of the 
computed value. These errors can be attributed to beam spillover, absolute calibra- 
tion errors, water surface inhomogeneities, and errors in the computed downwelling 
brightness field. 

5. Physical Basis for Polarimetric Signatures 


Insight into the relationship between various geophysical parameters and 
variations in the upwelling polarimetric brightness can be obtained by writing the 
upwelling surface brightness in terms of the bistatic scattering function [Tsan<7 
et ai, 1985]: 

T a {0 a ,<f> 3 ) = - — 4>*i fa) <i>i) sin 0id0id<j>i + T w [1 — n a (6 t , 4>,)} 
in p Jo Jo 

(28) 

where the surface reflectivity r a (6 t ,<f> t ) is: 

r Q {6.,4> t ) = / I ' 7/3 sin 0 t d6 t d<f>i (29) 

in p Jo Jo 

The downwelling background brightness can be considered to be the sum of an average 
component weighted over all incident angles and an angularly varying component: 


where 


= <m,«) + 


<r(*„*)> = J-E I’'' sinOMdh 

47r p Jo Jo 


(30) 


(31) 


With reference to (27), the second-harmonic azimuthal Fourier coefficient can now be 
written as: 


B a2 (0 t ) = 2 


[(r(0„<k)> - T w )r a2 (e.) 

2 * tv/2 


1 t 2 * tw/2 

+ — 5Z / / 4>i) 70a2(0 * j <f>i) sin 0id0id<j>i 

4n a Jo Jo 


(32) 


where 


r on (e.) = i~ r r 0 (e„ <f>.) d<j> t (33) 

Z7T JO 

is the n th azimuthal Fourier coefficient of the surface reflectivity, and 

-rpUe.AM - ^ jT d4,. (34) 
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is the n th azimuthal Fourier coefficient of jp a with respect to <f>,. 

The second harmonic amplitude B a 2 is influenced by two independent geo- 
physical effects, described by the two terms in (32). The first effect results from 
variations in the total surface reflectivity r Q (or equivalently, the surface emissivity) 
with fa, and is described by [(T(0i,fa)) - r u ,]r t , 2 (^). This effect causes surface 
emission and reflection harmonics in T a which are by nature of opposite phase. The 
second effect results from both elevational variations in the angular distribution of 
the downwelling Tp{&i, fa) and azimuthal variations in the bistatic scattering function 
Ipai and is described by the integral in (32). This effect causes background-induced 
harmonics in T a . 

Inspection of (32) shows that the amplitude of the second harmonic contains 
a component which is proportional to the difference between the weighted-average 
downwelling brightness and the water temperture. (This can easily be shown to be 
true for all harmonic orders.) Thus, a brightness contrast between the background 
and the surface temperature must be present to produce either surface or background 
harmonics. For example, if the background brightness was uniform over angle, un- 
polarized and identical to T w then the upwelling surface brightness T a would be 
unpolarized and equal to T w at all angles 4> t . 

The computed sensitivities of B o2 to various parameters provide additional 
insight into the link between geophysical parameters and radiometric observables. 
Accordingly, the numerical derivatives dB a7 /d(%SRH), dB a7 /d(h/ A) and dB a7 /dT w 
are shown in Table 1 for three values of 0, : 15°, 40° and 65°, representative of the three 
characteristic ranges discussed in Section 4. The calculations use quiescent surface 
and atmospheric states identical to those of Figs. 7a and b. 

Both positive and negative values of dB o2 /d(%SRH) are possible in clear 
air. This is explained by analyzing the components of B q2 under small increments 
in the background humidity profile. For 0, = 15°, 40° and 65°, the surface harmonic 
component in (32) is reduced due to an overall warming of the background brightness 
profile. In general the background warming also reduces 6T(0„ fa), however, this 
reduction is not uniform over the entire sky: elevation angles near the horizon are 
increased by a larger amount than those near zenith. This phenomona is characteristic 
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of the stratified nature of the atmosphere and varies somewhat with the particular 
atmospheric quiescient state. Since 't 0 a i{O g ,Oi,<t>i) exhibits both positive and negative 
values over 6 i, it is possible for the background-induced harmonic component in (32) 
to increase with a general warming of the background brightness profile. This is the 
case for all polarizations and for all three observation angles studied. 

For 0 g ~ 30°, the n = 2 surface emission harmonic dominates both the 
surface scattering and background-induced harmonics. Due to destructive interference 
of these harmonics, an increase in SRH decreases B a ?, as shown by the negative 
sensitivities in Table 1. When 6 t is increased to ~ 40° surface emission continues 
to dominate f? v2 , however, the surface scattering and background-induced harmonics 
dominate B ^ since the reflection coefficient is larger for a = h. Thus, the sensitivity 
values in Table 1 for a = h are positive. For 9 g = 65°, the n = 2 background- 
induced harmonics become large and dominate jB o2 for all polarizations, hence again 
the positive sensitivity values in Table 1. The background-induced harmonics are 
large when 0 S = 65° due to the rapidly changing background profile near the horizon 
along with the diffuse reflecting features of the sinusoidal surface (as described by 
iPai) when the observation angle is near the Brewster angle of the calm water surface 
(~ 71°). 

6. Implications for Remote Sensing 

The radiometric behavior under the GO model is a consequence of the slope 
distribution of a sinusoidal surface, the variance of which is — (y/2 irh/ A) 2 ss 0.049. 
Consider spaceborne passive remote sensing of ocean wave direction. For the open 
ocean, the slope distribution can be approximated using the model of Cox and Mxink 
[1954; Wilheit , 1979], which has the following variance: 

a\ = 0.003 + 0.0048 w (35) 

where w is the wind speed at a height of 20-m above the surface. Using w = 7 m/s ec 
(the most probable value) we obtain a 2 « 0.037, a value comparable to <r 2 . Thus, az- 
imuthal brightness signatures of comparable amplitude can be expected over striated 
portions of ocean. For 6 g outside of the range 35 — 45" the signatures are expected to 
be large enough for remote sensing of ocean wave direction, particularly for 0 g ~ 50°. 
However, signature amplitudes can be expected to be reduced somewhat by ocean 
foam, non-directional wave spectra and horizontal spectral inhomogeneity. 
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Since the azimuthal signature for Ty is in phase quadrature with both T v or 
T k% measurements of T v are expected to facilitate retrieval of ocean wave direction 
modulo 180° by reducing the ambiguity in the quadrant of the wave direction from 
fourfold to twofold ([e.g., Etkin et a/., 1993]. Since a real ocean surface is not sym- 
metric the Fourier amplitudes B an for odd n are nonzero. This was shown by Wentz 
[1992] for q = v, h using SSM/I data, and explained by an excess of leeward-side foam. 
Thus, elimination of the remaining 180° ambiguity is expected to be facilitated by 
the natural ocean wave asymmetry. In principle, ocean wave direction should be un- 
ambiguously measurable from single polarimetric observations, preferably at a large 
incident angle. In practice, the presence of foam, non-directional spectra and spec- 
tral inhomogeneity will likely necessitate additional observations at several azimuthal 

angles. 

Although the downwelling brightness field and the water vapor opacity be- 
tween the surface and a space-based observer (~ 3 dB) afTect the characteristics of 
the azimuthal radiometric features, signatures useful for ocean wave direction sensing 
at millimeter-wave frequencies (at least a few Kelvin in amplitude) should be ob- 
servable from space in clear air under nearly all humidity conditions and under light 
( < 0.1 kg/m 2 ) cloud cover. Indeed, the downwelling brightness depends on both the 
columnar water vapor content of the atmosphere and the amount of cloud water. At 
90 GHz and under US standard atmosphere conditions the clear-air zenith bright- 
ness ranges from ~ 15 K for 0% SRH to ~ 75 K for 100% SRH, and can increase to 
~ 145 K under saturated tropical summer conditions. In any of these cases the bright- 
ness contrast between the background and surface is large enough (~ 150 - 260 K) 
to produce observable signatures. For somewhat heavier clouds (0.1 g/m 3 from 1 km 
to 5 km altitude) the zenith brightness increases to ~ 190 - 220 K. Under such con- 
ditions the brightness contrast is reduced to ~ 50 - 80 K, and the amplitude of the 
azimuthal Fourier harmonics are reduced to ~ 30% of their clear-air values. In this 
case azimuthal brightness signatures might be observable from low-flying aircraft, but 
cloud opacity (~ 6 - 8 dB) would render them practically unobservable from space. 
The impact of clouds on spaceborne surface measurements is considerably reduced at 
lower frequencies (e.g., ~ 18 and ~ 37 GHz), albeit with reduced spatial resolution 
when using diffraction limited apertures of fixed size. 

As shown in Table 1 by the sensitivity to the parameter (h/ A) the Fourier 
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amplitudes generally increase with slope variance, and are particularly large for 6, 
exceeding ~ 50 — 60°. The anisotropic effects of striated water surfaces should be 
most readily observed at such angles. Otherwise, the sensitivity to T w is relatively 
small near 90 GHz. In addition, no significant changes occur upon incorporating salt 
at normal oceanic levels (3.5%). 

The results of this study are also applicable to nadiral- or near-nadiral sound- 
ing of atmospheric parameters. For angles 6, ~ 35°, the GO model suggests that water 
waves with slope variances comparable to open ocean can cause peak-to-peak varia- 
tions in the surface brightness of ~ 3 K, resulting in random variations of ~ 1 —2.5 K 
in a 90-GHz window channel. Comparable variations have been shewn to occur at 20- 
and 37-GHz over ocean [ Dzura et al ., 1992]. The magnitude of these variations are 
large enough to warrant consideration in algorithms for water vapor and temperature 
sounding, atmospheric wet-path delay measurements, and in satellite climatological 
studies of parameters that might be statistically correlated with ocean wave or surface 
wind direction. 

7. Summary 

Presented in this study are the results of controlled partially polarimetric 
measurements of thermal emission at 91.65 GHz from a striated water surface as cor- 
roborated by a geometrical optics radiative model. The measurements were obtained 
outdoors using a precision polarimetric radiometer which directly measured the first 
three modified Stokes’ parameters. Significant variations in these parameters as a 
function of azimuthal water wave angle were found, with peak-to-peak variations in 
Ty of up to ~ 10 K. The measurements are well corroborated by the GO model 
over a range of observations angles from near nadir up to ~ 65° from nadir. The 
model incorporates both multiple scattering and a realistic downwelling background 
brightness field. 

Both the data and the GO model suggests three characteristic ranges of 
observation angle. The largest azimuthal signatures in T v , 7\ and Ty are found at 
observation angles beyond ~ 50° from nadir. For 35° ~ 0, ~ 45° the azimuthal 
signatures practically vanish. For 6 t < 35° the azimuthal signatures are of smaller 
amplitude, with phase reversals occurring in T v and Ty. The presence of water waves 
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was also noted to influence the azimuthally-averaged values of T v and 7^ by several 
Kelvins. Applications of the study include passive polarimetric remote sensing of 
ocean wave direction from space, oceanic wave studies from both space and aircraft, 
and prediction of polarimetric surface effects for atmospheric sounding. 
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Table 1: Computed sensitivities dB a2 /dx for x = %SRH,{hj A) and T w , and 
v,h,U. 





dB a2 jdx 


Dimen- 

X 

e. 

a = v 

a = h 

Q=U 

sions 


15° 

-0.014 

-0.016 

-0.030 


%SRH 

o 

O 

-0.030 

0.016 

-0.052 

(K/%) 


O' 

o 

0.021 

0.168 

-0.056 



15° 

63.2 

53.3 

116.9 


h/A 

40° 

22.1 

82.4 

-29.1 

(K) 

65° 

213.2 

304.2 

319.7 



15° 

0.009 

0.006 

0.014 


T w 

40° 

0.018 

0.002 

0.016 

(unitless) 


65° 

0.054 

0.045 

-0.003 




22 


Figure 1: Block diagram of the 3-channel (v*, vg and vc ) polarization-correlating 
radiometer. 


Figure 2: Experimental configuration for polarimetric radiometric measurements of 
emission from a striated water surface. 


Figure 3: Geometry associated with the GO brightness model. 


Figure 4: Multiple scattering and shadowing thresholds for a sinusoidal surface z = 
h sin( 27 rx/A). 


Figure 5: Geometry associated with multiple scattering from a sinusoidal surface. 


Figure 6: Computed upwelling brightness perturbations for 0° < 9, — 68°, 

0° < <f> t < 90°, hj A = 0.05 and a surface relative humidity of 50%: (a) T v — Tj; c \ 

(b) T h - Ti‘\ (c) T V . 


Figure 7: Azimuthal Fourier coefficients for the brightness perturbations in Figs. 6a-c: 
(a) DC Fourier amplitude, (b) 2 nd harmonic magnitude i? o2 , and (c) 2 nd harmonic 
phase $ a2 . 


Figure 8: Measured and calculated polarimetric brightnesses for 9, = 20°, T w = 
291 K, h/A = 0.05 and SRH = 85%: (a) T v and T h , (b) T v . 


Figure 9: Measured and calculated polarimetric brightnesses for 6, = 39°, T w = 
291 K, h/A = 0.05 and SRH = 80%: (a) T v and 7*, (b) T\j. 


Figure 10: Measured and calculated polarimetric brightnesses for 6, = 65°, T w = 
288 K, h/A = 0.05 and SRH = 50%: (a) T v and T h , (b) T v . The dashed line in 
(a) and (b) shows the computed values of Th and T\j (respectively) when multiple 
scattering is neglected. 
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